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ABSTRACT 


Advanced S.T,O.L. aircraft are critically sensitive to loss of 
engine thrust during landing, with the problem of encountering unaccept- 
ably high sink rates. The problem is aggravated by the high degree 
of stability augmentation common to this class of aircraft which denies 
the pilot the immediate motion cues required to manually arrest large 
transient sink rates. Automatic control is required to arrest the 
immediate transients, allowing the pilot time to recognize the situa- 
tion and decide upon appropriate action. 

Modern optimal control theory, which determines state variable 
feedback controllers for multivariable systems, is especially tailored 
for application to this S.T.O.L. control problem. System performance 
can be specified in terms of the desired response characterized in the 
system state variables. The resulting performance is typically better 
than that obtained from classical design procedures. 

The purpose of this research is to define candidate autopilot 
control laws that control the engine failure transient sink rates, by 
demonstrating the engineering application of modern state variable con- 
trol theory. 

This work provides a comparison of the results of approximate modal 
analysis to those derived from full state analyses provided from com- 
puter design solutions. The aircraft is described and a state variable 
model of its longitudinal dynamic motion due to engine and control 
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variations is defined. The classical fast and slow inodes are assumed 
to be sufficiently different to define reduced order approximations of 
the aircraft motion amendable to hand analysis control definition methods. 
The original state equations of motion are also applied to a large scale 
state variable control design program, in particular OPTSYS. The re- 
sulting control laws are compared with respect to their relative res- 
ponses, ease of application, and meeting the desired performance ob- 
jectives. The limitations of minimizing feedback paths is investigated 
with the objective of utilizing those currently in use by the classical 
autopilot control law. The performance of each resulting control law 
is demonstrated by digital computer simulation. The sink rate transients 
are shown to be controlled within 5-10 sec after an engine failure and 
the peak sink rates are shown to be less than 19 ft/sec. 
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I. INTRODUCTION 


1.1 STATEMENT OF ^ PROBLEM 

The critical phases of aircraft flight are landing and takeoff. 
This is especially true of short distance takeoff and landing (S.T.O.L.) 
transports because they fly slowly during these phases depending upon 
powered lift augmentation to maintain the desired flight path. In this 
thesis, the landing phase is studied of S.T.O.L. flight of aircraft 
using lift augmentation from the primary engines, typified by the YC-14 
and YC-15 transports currently being flight tested. It is the purpose 
of this study to minimize the effects of loss of one engine during the 
final phases of approach (less than 500 feet altitude). The problem 
has been defined in the similar research of Messrs. Franklin and 
Nieuwnhuisse at NASA-Ames (N-1), wherein several simple solutions were 
attempted without the desired results. 

The approach to this problem is to utilize the integrated analysis 
feature of modern state variable control theory to define several can- 
didate feedback control laws that minimize the buildup of sink rate when 
an engine fails. The relative perfoirmance of these control laws is 
compared by linear simulation. The effectiveness of the control laws 
is demonstrated. 

Preliminary familiarization efforts with the study aircraft and 
engine performance bri-ngs about the conclusion that the engine response 
is significantly faster than the expected closed loop airframe response 
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that is required to control major sink rate disturbances. Therefore, 
the approach used in development of a workable control law is to per- 
form design analyses without considering the engine response constraints. 
Engine influence is evaluated with the control law to determine the 
validity of this assumption by simulation studies. Two control analyses 
methods are researched to identify the relationship between standard 
aircraft modes and the dynamic modes of the study aircraft configura- 
tion. A byproduct of this approach is the comparison of approximate 
analyses methods (valid when significant mode separation exists; such 
as is the case in convention aircraft) to optimal control algorithms 
which account for full state interactions. Second order effects as well 
as resolving resulting engine influence problems were not included. 

Topics for further research associated with these issues have been 
identified, 

1.2 ENGINE FAILURE DURING APPROACH 

When an engine fails during the landing approach of a S.T.O.L. 
aircraft, excessive sink rates may occur before the pilot becomes aware 
of the failure. 

The aircraft of interest in this study obtains lift augmentation 
through directing the engine flow over the extended flaps. Loss of one 
engine immediately reduces the lift by 10% with an increase in drag due 
to the engine nacelle and the asymetric trimmed flight condition. The 
aircraft is control augmented with pitch and roll/yaw stability aug- 
mentation and a yaw trim control. These systems respond so as to main- 
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tain the aircraft heading and minimize transients. The aircraft drag 
is increased due to the control deflection and asyraetric flight 
attitude required to maintain the proper heading. 

All of these actions, on the part of the automatic controls, 
virtually eliminate any immediate cues to the pilot that an engine has 
failed. The only remaining cue is the resultant sink rate due to loss 
of an engine without a compensating increase in power setting. Separate 
engine failure indicators provided to the pilot (researched in Reference 
N-2) do not alleviate the sink rate because the pilot's workload keeps 
him too busy to notice the indications, especially during turbulent 
approaches. Therefore, an additional automatic control loop is required, 
either internal to the S.A.S. or specially triggered by an engine fail- 
ure. 

1.3 CONTROL STUDIES OVERVIEW 

The control studies are directed to determine state feedback con- 
trol laws that arrest the sink rate due to the loss of one engine. The 
feedback control is defined using the five state plant model as a basis 
for development. Once a particular control law is defined, it is tested 
in simulation with the engine model adjoined to the plant model. This 
approach is employed because the engine performance model provided for 
this study responds slightly faster than the fast airframe dynamics 
(this is expanded in a subsequent section). Two analytical approaches 
are employed to contrast the results and gain understanding of the 
"roles" of each of the controls in relation to the states. 
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Three control concepts are investigated to determine their relative 
effectiveness in controlling sink rate and maintaining acceptable air- 
craft response. The bulk of the work utilizes the three controls (which 
seem to be sufficient in the proper mixture) of thrust, elevator and 
aft flap. The spoilers, nominally deflected 30 degrees, are added to 
the control list to determine if fast direct lift control significantly 
increases the response in arresting sink rate. These control options are 
employed using full state feedback which includes the altitude state. 

The third concept studied removes the altitude state feedback to deter- 
mine if sufficient damping is available to eliminate excessive sink 
rates. 

The first control concept is analyzed using reduced order approxi- 
mations to define feedback control laws. The standard approach of 
separating the short period and phugoid/altitude modal equations of 
motion is attempted. Such methods work well if these modes are widely 
separated. The separation of modes in this case is marginal making it 
an interesting attempt to obtain workable solutions where the, poles are 
situated so they may interact severely. Two solutions are generated 
from this approach. 

The value of the approximate study is not wholly in determining 
working control laws. The analysis provides experience with the relation- 
shipment of airframe response eigenvalues and the controls. The Root Square 
Locus method (Reference B-3) of defining optimal approximate solutions 
provides a format to observe the relative power of controls in providing 
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desired responses and also observe the root interactions that may cause 
trouble. 

The second analytical approach employs the full five state plant 
model and the power of automated optimal solutions provided by the Stan- 
ford OPTSYS Computer Program (Reference B-1). Each control concept is 
analyzed with this method, obtaining control law results for simulation. 
This method has one drawback, however, it assumes full state feedback, 
since it is common that some states are not available for feedback. 

The standard solution to this is to define an estimator to provide the 
necessary states. Another approach is to eliminate the unmeasured state 
feedback from the control law without recomputing the remaining feedback 
gains. Since this latter scheme is economical it is investigated as 
part of the control analysis. 

1.4 ENGINE MODEL ANALYSES 

Engine throttle response curves of the study aircraft's power 
plant are provided. These are studied to define a first order response 
model for simulation analyses. The previous work included three charac- 
teristic engine models that had thrust transient responses ranging from 
slow to fast. (Reference N-2). The preliminary results from the engine 
performance provided herein show this engine to be about equally res- 
ponsive as the medium response engine model from that previous study. 
(Reference N-2). 

A closer look using identification analyses to define higher order 
models to determine their effects (slower roots perhaps) is identified 
as a suggested topic for continuation of this work. 


1.5 SIMULATION MODEL WORK 

A linear simulation model is provided in FORTRAN IV to facilitate 
demonstration of the control law results and comparison of state devia- 
tions and control deflections resulting from the control laws. Several 
of the control law results are simulated and their relative performance 
is compared, 

1.6 CONTROL LAW APPLICATION 

The feedback control laws developed in this study are perturbation 
control laws, whereas the full state values are measured. The con- 
version of these control laws into an aircraft application is discussed, 
with specific attention to the need of unmeasured state feedback. The 
study aircraft is defined to have speed, pitch, and pitch rate readily 
measured. Altitude is also available as either a direct measure or an 
open loop observer that is said to work as well as a direct measure. 
Angle of attack is not measured. 
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II. COORDINATE SYSTEM 


The coordinate system employed is a right handed orthogonal basis 
presented in Etkin and Seckel (References E-1 and S-1) modified to pro- 
vide positive parameters defining sink rate. The positive pitch angle, 
theta (0) is down and sink rate h is also positive down, allowing the 
sink rate equation to be made up of a sum of positive quantities. This 

f 

requires theta rate (©) to be opposite in sign to q (the latter being 
defined in the standard convention of positive nose up) . This coor- 
dinate system is shown on Figure II-l indicating positive quantities of 
the states and forces on the aircraft. The intertial and body fixed 
axes are sho\m with the pitch axis directed out of the page. The sys- 
tem equations of motion are derived in these coordinates in Appendix A. 
Where results are shown (such as the simulation state histories in 
Chapter X, the axes of the modified states are inverted making the 
visual association of downv/ards orientations being below the axes origin 
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III. AIRCRAFT CONFIGURATION 


3.1 SUMMARY 

The particular S.T.O.L, aircraft chosen for this study is a four 
engine high wing transport which uses an externally blown flap system 
for lift augmentation. The characteristics of the study aircraft are 
described in Reference N-2 and N-3. The particxilar details pertaining 
to this research (longitudinal dynamics) are repeated here. The air- 
craft is a study model currently in use at NASA Ames in moving base 
simulator studies. 

3. 2 STUDY AIRCRAFT DESCRIPTION 

The study aircraft configuration is displayed on Figure III-l, 
and its longitudinal physical characteristics are summarized on Table 
,III-1. It is a 150,000 lb transport similar to the McDonnell-Douglas 
YC-14 currently being flight tested. The lift augmentation is provided 
during takeoff and landing approach by externally blown flaps, as sho’t'm 
on Figure lll-l. Unlike conventional aircraft, the flaps extend into 
the engine efflux. The high speed flow is diverted by the flap struc- 
ture, increasing the circulation over the wing and increasing the lift 
provided by the v;ing. 

The nominal approach configuration is used in this research. The 
nominal approach is trimmed for equilibrium flight with the flaps exten 
ded to the 60 deg position. 
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The study aircraft in the approach configuration has four 
longitudinal controls available. These longitudinal controls are: 
elevator, engine thrust, spoilers, and aft flap. The elevator operates 
the same as conventional aircraft; the elevator controls pitch attitude 
and rate. (The thrust or throttle modulates lift and drag.) The 
spoilers are usually operated differentially for lateral control. 

However, with the spoilers initially deployed and operated together 
they can provide direct lift control, with response faster than the 
typical engine throttle response. The aft flap control consists of 
the last element of the flap structure, which is geared to move rapid- 
ly like the elevator or spoiler surfaces (sho^-m on Figure lll-l). Where- 
as the main portion of the flap is slow moving and considered to be a 
variable configuration device deployed to varying degrees when low speed 
and high lift is required. The aft flap is principally a drag modula- 
tion control. 
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TABLE III-l 


S.T.O.L. TRANSPORT PHYSICAL CHARACTERISTICS 

Weight 150,000 lb 

Pitch Moment of Inertia (lyy) ..... 2,105,000 slug ft 

Wing Area (AREF) 1,667 ft^ 

Mean Aerodynamic Chord (CREF) ..... 16.3 ft 

Horizontal Tail Area 570 ft^ 

Maximum Engine Thrust 84,000 lb 

(Sea Level Static) 



IV. AIRCRAFT AERODYNAMICS 


4.1 SOURCE OF AERODYNAMIC DATA 

The aerodynamic coefficients for this study are provided in raw 
form. The data, given in Reference N-3, requires modification before 
computing the numerical derivatives in the equations of motion. Most 
of the longitudinal coefficients are given in standard form referenced 
to the wing area (AREF) and mean aerodynamic chord (M.A.C.). However, 
two coefficients, the lift and drag due to elevator deflection are 
given referenced to the horizontal tail area in Reference N-3. 

4.2 AERODYNAMIC DATA MODIFICATION 

The raw form aerodynamic data must be modified for use in this 
study. The raw data is not in a consistent form, two aerodynamic coef 
ficients must be modified to be consistent with the reference area 
(AREF). The data is digitized for use in this study to allow computer 
solution of the dynamics matrices in state variable form. 


The coefficients 
referenced area: 


Ct and Cr, must be normalized to the V7ing 
de de 


Cl = Cl* 

^de ^de 


Se 


^Tail 

AREF 

^Tail 

AREF 
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The aerodynamic coefficients of powered lift aircraft include the 
thrust effects on lift, drag, and moments, and it is nearly so in this 
case. The ram drag component due to the air entering the engine nacelles 
is not included in the coefficients of drag (CD) and must be added se- 
parately. Ram drag is part of the forward thrust contribution of the 
engines, where the thrust is due to the air mass accelerated through 
the engine (Reference D-1). 

T - ('^exhaust “ ^aircraft^ 

T - ^ Q Vaircraft 

T - Vexhaust " Uq 

The latter quantity is ram drag where is mass flow rate. 

Dividing by reference area (AREF) and dynamic pressure (q) provides the 
drag coefficient contribution due to ram drag. 

^CD gQ ’Jq 
q AREF 

4.3 STUDY AERODYNAMIC COEFFICIENTS 

The raw data for the 60 deg flap deflection configuration (given 
in Reference N-3) is digitized into a consistent form referenced to the 
wing geometry. These data are shown in tabular form in Appendix C 
from the data in Reference N-3. The following coefficients are pro- 
vided in graphic form as a function of angle of attach (<^) and thrust 
coefficient (cj) in Reference N-3, 
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CL 

CD 

CM 

''I^de 

Cn 

®de 

Se 

%f 

Cn 

®df 



The apparent plotted data points are read to provide the digitized data, 
to minimize the introduction of significant interpolation errors. However, 
the spoiler data is provided in table form. The coefficients CL^jgp and 
CD(jsp are given over a range different from the plotted data. CM^sp 
not given, so it is assumed the spoiler deflection does not affect the 
aircraft moments. These table data are interpolated to provide digital 
points at the same conditions the plotted data are taken for presentation 
in appendix c and use in this study. 

The aerodynamics are not modified for ground effects. The study 
is done without ground effects, leaving their contribution for evalua- 
tion by simulation, as the resulting control must work out of ground 
effect as well. Ground effects increase the lift near touchdown so 
their elimination is a conservative choice. 
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V. EQUATIONS OF MOTION 


5.1 FOBM OF THE EQUATIONS OF MOTION 

The plant dynamic equations are derived in state variable form, 
in Appendix A as linear dynamic equations about unperturbed (equilibrium) 
steady approach condition. The resulting linear state variable model is 
of the form: 

X ~ Fjc + Gii 


In particular the equations are; 


• 

u 


"Xu 

^w 


- w 

q 

o“ 


u 


Xcj 

de 

^df 

^dsp 


dcj 

w 


Zu 

Zw 


- u 

0 

0 


w 


Zcj 

Zde 

Zdf 

Zdsp 


de 

q 


Mu 


^^q 
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0 


q 


Mcj 

Mde 

Mdf 

Mdsp 


df 

e 


0 

0 

-1 


0 

0 
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0 

0 

0 


dsp 

ji 


_0 

1 

0 


u 

0 


h 


0 

0 

0 

0 




The derivatives are defined in terms of the nominal states and 
the airframe aerodynamic coefficients in Appendix B. 

The plant model consists of the equations of motion with the dimen- 
sional derivatives computed about the nominal approach condition for the 
particular configuration of the plant model requires the equilibrium 
trimmed approach conditions be first determined, from which, the dimen- 
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sional derivatives are computed. The approach employed is to generate 
an equilibrium performance map for each configuration of interest;, 
select the desired approach condition; and compute the dimensional 
derivatives. 

Two computer programs are built to facilitate computation of the 
derivatives, one which computes a map of equilibrium flight conditions 
for a specified configuration, and another that computes a particular 
equilibrium flight condition and the corresponding perturbation de- 
rivatives in the above dynamic equations. These assist the designer 
in providing numerical forms of the equations of motion, directly in 
the form for analysis and/or simulation, for any equilibrium approach 
condition and configuration of interest. The latter of these two pro- 
grams is used in conjunction with a linear simulation model to demon- 
strate the resulting control performance, providing the capability to 
evaluate control laws at off design conditions by specifying the 
appropriate equilibrium conditions without recomputing the dynamics 
matrices. 

5.2 EQUILIBRIUM PERFORMANCE MAP 

The equilibrium performance map is derived for the range of angles 
of attack from -8 deg to 20 deg, and over the thrust range of the study 
aircraft. The configuration is defined by selecting the nominal position 
of the aft flap and spoiler controls (with the main flap deflection 
60 deg). The force equations are solved by iteration using the digitiz- 
ed aerodynamics tables. The particular equations are: 
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CM = 0 (eg trimmed condition) 

q = % 

cj = T/(q.AREF) 

T = TAN-1 ( -CD/CL) 
g*m*COS()f) = CL*q*AREF 

Where the aerodynamic coefficients represent the sum of aerodynamic 
derivatives in the form: 

CL = CL (o<:, cj) -y CL^g (o<,cj) de + CL^^ ( c<, cj) df + CL^gp ( «< , cj) dsp 

An iterative algorithm is programmed in FORTRAN IV to provide the 
desired performance maps. This program is included in Appendix D. 

Equilibrium performance maps are generated for the two configura- 
tions considered in this study to indicate the magnitude of the problem 
of engine failure during approach. The nominal approach configuration 
is with the flaps deflected 60 deg and the spoilers retracted. The 
nominal approach condition of approximately 80 kt provides a 5*6 deg 
flight path angle with the elevators trimmed at -3.2 deg deflection 
trimming the aircraft at 2.2 deg angle of attack (Figure V-1). The 
nominal approach sink rate is 13.2 ft/sec. The nominal approach thrust 
is 50400 lb at 60 percent of the maximum thrust of 84000 lb. A loss 
of one engine provides a thrust of 37800 lb. The equilibrium flight 
condition at 80 kt with one engine out is trimmed at 6 deg angle of 
attack providing a flight path angle of 8 deg (a sink rate of 19 ft/sec. 
If the initial trim angle of attack of 2 deg is maintained, the flight 
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conditions are 87 kt at a flight path of 7 deg (an 18 ft/sec sink rate). 
The constant thrust curves show the aircraft to be operating on the back . 
side of the power curve where increased angle of attack results in in- 
creased sink rate. It is also interesting to notice that the full 
throttle capability of this aircraft in the approach configuration with 
one engine failed (63000 lb of thrust) does not provide level flight, 
but does provide a significant margin to maintain the approach conditions 
(Figure V-1). Therefore, the configuration must change, by retracting 
the flaps, to perform a go-around (e.g. abort the landing). 

The possibility of approaching at a slightly higher speed with the 
spoilers deflected is considered in this study. A 6 deg approach is 
possible at 84 kt with the spoilers deflected 30 deg at the nominal 
approach thrust and angle of attack (Figure V-2) . In this condition, 
it is expected that the control response upon engine failure be to re- 
tract the flaps. Viewing the equilibrium conditions with flaps retract- 
ed at 84 kt and one engine failed shows the flight path to be 7.5 deg 
which has a sink rate of 18.5 ft/sec (Figure V-1) at two degrees angle 
of attack (requiring an additional 4 percent of the maximum total thrust 
to maintain the flight). The flight path is 6.5 degrees and a sink 
rate of 16 fps. The spoilers may alleviate the sink rate somewhat but 
additional control and throttle response is required to maintain a 
nominal approach flight path. 
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Y (deg) 


CONFIGURATION 

4f = 60 deg 
df = 0 deg 
dsp - 0 deg 



T (deg) 


CONFIGURATION 

= 60 deg 
df = 0 deg 

dsp = 30 deg 



5.3 COMPUTATION OF DIMENSION DERIVATIVES 


The dimensional derivatives that comprise the F and G matrices in 
the equations of motion are computed about an equilibrium approach con- 
dition. The method is to first determine the angle of attack of equi- 
librium for specified thrust, control settings, approach speed and 
flight path which are approximately detemined from the performance 
maps. The equilibrium solution here serves to determine the angle of 
attack condition to more accuracy than is apparent from visual inter- 
polation of the equilibrium approach performance curves, when a part- 
icular velocity and flight path is desired for study. 

Once the appropriate trim equilibrium conditions are determined, 
the dimensional derivatives are computed from the aerodynamic coef- 
ficients and the equilibrium flight conditions . The computational de- 
tails for the dimensional derivatives in the coordinates of this study 
in Appendix B and the FORTRAN V computer program is presented in 
Appendix D. 

The dynamics matrices are computed for the two configurations 
described in the previous section (Section 5.2). The nominal approach 
conditions are: (Where units are consistently in ft, rad, sec with the 

control deflections in deg) 

F MATRIX 


-.0372 

.0766 

-5.2114 

32.2000 

.0000 

-.4769 

-.4424 

135.2396 

.0000 

.0000 

-.0000 

- . 0046 

-.7457 

.0000 

.0000 

.0000 

.0000 

-1.0000 

.0000 

.0000 

.0000 

1.0000 

.0000 

135.2396 

.0000 

.0000 

.0000 

.0000 

.0000 

.0000 
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G MATRIX 


.9102 

-.0269 

-.2230 

.0002 

9.6367 

-.0899 

-.2631 

.0547 

-.0393 

-.0261 

.0133 

.0000 

.0000 

.0000 

.0000 

.0000 

.0000 

.0000 

.0000 

.0000 


Which are of the form described in Section 5.2, 

The dynamics matrices at the approach condition with 30 deg 
of spoiler deflection are; 

F MATRIX 


0393 

.0730 

-5.2933 

32.2000 

.0000 

4564 

-.4222 

140.4003 

.0000 

.0000 

0000 

-.0051 

-.7647 

.0000 

,0000 

0000 

.0000 

-1.0000 

.0000 

.0000 

,0000 

1.0000 

.0000 

140.4003 

.0000 


G MATRIX 


.9652 

-.0279 

-.2295 

.0001 

10.9655 

-.0952 

-.2678 

.0579 

- . 0422 

-.0279 

.0146 

.0000 

.0000 

.0000 

.0000 

.0000 

.0000 

.0000 

.0000 

. 0000 _ 


The two nominal approach conditions have nearly the same dynamics matrices. 
This similarity provides the basis for an expediate assumption applied to 
the spoiler analyses (discussed later in this report) , 
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VI. APPROXIMATE MODEL ANALYSIS 


6.1 THE METHOD OF MODAL ANALYSIS 

The standard aircraft approach of separating the analysis of 
short period and phugoid modes is investigated. This results in re- 
duced order equations that are amendable to hand analysis techniques 
such as successive loop closure. The practice of analyzing the two 
basic aircraft modes independently is based upon spectral decomposition 
concepts. When the characteristic foots are widely separated, then 
their characteristic eigenvectors, which are composed of combinations 
of the states that are disturbed by each mode (set of characteristic 
roots), exclude the states of other modes. Therefore, any particular 
mode is composed of a small subset of the aircraft's states, allow- 
ing analysis of that mode with the other modes undisturbed (i.e. 
perturbations remain zero) . The open loop characteristic roots are 
moderately separated in the case of the study aircraft (Figure VI- 1). 

If feedback control can move these roots further into the left half 
plane, while maintaining separation of the modal state responses, 
the approximate solutions will be valid. 

6.2 SHORT PERIOD MODE ' 

The approximate analysis first proceeds in the same manner as 
for conventional aircraft; beginning with the short period (fast) 
mode to obtain increased damping; then approximate phugoid- altitude 
mode equations are defined to develop stable responses. The short 
period (made up of w and q) has very little feedback of other states 
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s- PLANE ( 

1 - 2 , 



FIGURE VI-1 AIRCRAFT OPEN LOOP POLES 


(u, 0, and h) and feeds those states derivatives. The partioned short 
period equations are; 


*■ -s 

w 

q 

u J 


-.4420 135.2 

-.0046 -.7461 


-T r 

w 


-.090 

-.026 


de 


choosing the elevator control for moment damping. 


Assume a feedback control of the form 

X = Fx + Gde 
de = (0, Cq) x 

to increase the system damping. The resulting characteristic equation 
is determined from: 


F - GC = -.442 135.2 + .09 Cq 1 

‘ I 

-.0046 -.7461 + .026 CqJ 


which provides the characteristic equation from the .determinant of the 
dynamics matrix; 


s - (F - GC) = 0. 

Resulting in the characteristic equation; 

3 ^ + (1.188 -.026 Cq) s + .9517 -.0111 Cq = 0 

The root locus for this system is sketched as a function of Cq on Figure 
VI-2. A gain of -47.4 deg sec/rad provides roots on the real 
axis. The study by Franklin shows that the existing gain selection is 
-40 deg sec/rad which provides poles at s = 1.116 + j .389 sec"^. This 
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latter gain is selected since it provides acceptable damping and is cur- 
rently in use in other studies of this aircraft. The approximate phugoid 
equations are defined with the short period feedback loop closed, 

6.3 PHUGOID AND ALTITUDE CONTROL 

The approximate phugoid equations are determined by making the 
quasi-steady assumption that w and q are zero. Then, solving for w and 
q in terms of the other states (u, 0 , and h) and the controls (dc j , de, 
and df), so that w and q can be replaced in the other state equations. 

The relationship between vr and q to the other states is determined from 
the w and 4 of the full state dynamical equations: 



/-0.4772 

I 0 


\!ul + 

^0.442 131.6 ' 

w 

•f /-9.64i -0.2532 0.05691 

dcj 

i J 

■ 


1 

df 

/ 

-0.0046 -1.1790, 

q 

, » 

[-0.0391 0.0133 0 / 

dsp 


which results in the following relationship for w and q: 


• * 

w 

= 

'-0.6117 ^ 


q 


0.0015j 



-16.04 

0.0194 


0.9159 0.073 

0.0051 ■ -0.0002 


dcj 

df 

dsp 


replacing w and q in the equations for u,0 , and h results in the approxi- 
mate phugoid equations with the short period loop closed: 


r* 1 


u 


• 


e 

=: 

h 


- . 



-0.0919 


-0.0047 

-0.6117 


32.20 0 

0 0 


135.2 Oj 
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-0.4233 -0.1793 0.0072 


dcj 

+ 

-0.0194 -0.0051 0.0002 


df 


_-16.04 0.9159 0.073_ 


dsp 


These approximate equations are used with the two remaining controls 
dcj and df (dsp is investigated later). 

The remaining states u, 0, and h are to be controlled by dcj and 
df. The two controls allow only tv?o states to be directly addressed, 
depending upon the resulting system interactions to provide acceptable 
response of the remaining state. The method applied to determinine con- 
trol laws for dcj and df is successive application of the root square 
locus technique, which is described in Reference H-1. 

The speed and altitude ‘ states are chosen for control by df and dcj 
respectively. The 0 state is chosen to be implicitly controlled because 
it has the smallest coefficients in the dynamical equation and it is a 
state that contributes in a small way to the problem at hand. The choice 
of df to control u is due to the large impact of dcj on ^ ^i/ ddcj is 
38 times ^u/ ddcj while ^h/ idf is only 5 times more than ^u/ ^df). At 
this point it appears natural to connect dcj with h. 

The speed control loop is determined from the u/df transfer function 
in the s plane, which is: 

u_(s) = -0.1793 s (s -1 0.9159) 
df s (s2 + 0.0919 s -!- 0.1513) 
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The root square locus is formed about the optimal weight ratio 

u (s) Ay u (-s) +1 =0 

df df 

The locus is sketched on Figure VI-3. The closed boxes on Figure VI-3 
indicate the location selected for closed loop roots which correspond to 
A^ of 14.25 deg2 sec^Vft^. The roots are located at s - -.6 + j .525 

sec“lj to provide approximately 0.6 - 0.7 damping ratio. The feedback 
gains are determined from an assumed feedback configuration by matching 
the like coefficients of the transfer function and solving for the un- 
determined feedback coefficients. In this case two feedback gains are 
sufficient to' define the quadratic root location. A feedback law of 
the form: 

X = (F + GC) X 

is assumed, where the feedback coefficient matrix is assumed to be: 

C = (Cu, C@ , 0) . 

Entering this form of feedback into the approximate phugoid equations 
provides the following characteristic equation: 

s^ -h s2 -I- (0.0919 + 0.0051 Cq -f- 0.1793 Cu) s 

+ 0.15134 - 0.00037 Cq +0.1642 Cu = 0. 

which is to look like the following as determined from the root locus 
diagram: 

s (s^ + 1.2 s + .636) = 0. 
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FIGUPOi: VI-3 ROOT SQUARE LOCUS u (s) 



Solving for the feedback gains results in: 

C = (Cu, Cq , 0) = (3.19, 105.26, 0.) 

in consistent units of deg per ft, sec and rad. The equations for the 
next control detemination are: 

X = (F + GC) X + Gu 


where the F + GC matrix includes the feedback computed above. The 
analysis is repeated for the next successive control, in this case dcj 
which is chosen to control h. The nev/ dynamical equations are: 


• 

u 


-0.663 

13.327 

0^ 


u 


-0.423 ” 



-3.196 

0.5368 

0 


© 


-0.0194 

h 


2.306 

231.646 

0 


h 


-16.04 

. — 



Again, determine the transfer function: 

h (s) = - 16.04 (s^ i- .4676 X 4- 22.889) 

dcj s (s^ + 1 . 2 XT- . 636) 


which is put into root square locus form: 


h ^(s) ^h h (-s) -f 1 

dcj Bcj dcj 



The resulting root locus is sketched on Figure VI-4. The locus is 
sketched as a zero degree locus which provides a stable path along the 
real axis for the root at the origin (the h root) . Notice that the path 
of the real root is unbounded. However, this approximate analysis is 
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FIGURE vr-4 ROOT SQUARE LOCUS h (s) 

dcj 
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predicated upon the separation of modes and the short period mode roots 

lie about the s = -1.0 sec"^ region. The pole selection of the previous 

part (u/df) is half way between the origin and the short period territory 

and. possibly interacting. The selection of root locations here, is 

made to just make h stable and not move the other roots to far out. 

Additionally, the asymptotes of the oscillatory roots are at less damping 

than the initial root location, so the less they're moved the better the 

damping of the phugoid, barring major interactions. Selecting the real 

root location to be s = -0.2 sec“l meets the criteria. This location is 

obtained with a weight ratio of ^h = 0.0003 (ft“2), v?here upon the 

Bcj 

oscillatory roots move imperceptibly. The approximate characteristic 
equation is : 

(s + .2) (s2 -t 0.467 X + 22.89) = 0. 

The cubic requires three feedback coefficients of the form: 

C = (Cu, Ce » Oh) . 

Proceding as in the previous control case (equating coefficients of the 
closed loop characteristic equations) provides the following control 
gains defining dcj: 

C = (0.175, 13.762, 0.0125). 
in consistent units of ft, sec and rad. 

6.4 EVALUATION OT Tl^ CONTROL LAW 

The approximate modal analysis, just described, provides a set 
of control lav7S, obtained by ignoring any interaction between the tra- 
ditional aircraft modes. The resulting expected root locations are 
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indicated by X’s on Figure VI-5. The dotted lines indicate the expected 
travel of the unaugmented poles to the closed loop positions. The phugoid 
poles are moved out into the region occupied by the short period poles 
before augmentation and also close to their closed loop location. It is 
possible that significant root interaction has taken place, changing 
these results. This is tested by computing the closed loop poles of the 
full five state model with the control law just determined (with the aid 
of the OPTSYS Program). 

The closed loop pole locations of the five state model demonstrate 
that significant mode interaction does occur. The open boxes shorn on 
Figure VI-5 are the actual closed loop pole locations provided by the 
control law. 

The true results compared with those of the approximate sol- 
ution appear very different. The true poles are much slovrer than ex- 
pected and -the real pole is much faster; their relative real locations 
are switched about. The approximation assumption (which assumes no 
interaction) is violated severely; there being great apparent interaction 
between the w, q, u, and h states. The poor association of results 
provided by the approximate mode analysis occurs because the original 
slow mode is driven into the region of the fast modes. Therefore, if 
this fast and slov^ separation can be maintained, the validity of the 
approximation shall be also. This hypothesis is tested by arbitrarily 
locating the phugoid roots closer to the origin. 
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FIGURE VI-5 CLOSED LOOP POLES 

APPROXIMATE OPTIMAL SOLUTION 
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6.5 REVISED APPROXIMATE PHUGOID ANALYSIS 

The approximate analysis is revised to maintain the validity of 
the spectral mode separation assumption. Stepping back to the first loop 
closure (u/df) of the previous analysis, it is observed that the two os- 
cillatory poles can be placed arbitrarily with the two feedback coef- 
ficients C^ and . The trouble is caused by the oscillatory poles being 
moved out two far to the left, into short period territory on the s plane. 
The objective here is to increase the damping of these two poles without 
moving them into short period territory, as did the root square locus 
technique. The following choice of feedback gains meets the criterion 
by placing the oscillatory poles at s = -0,39 + j 0,176 (sec~l): 

= 0.4946 (deg sec/ft) 

Ce = 119.9 (deg/rad). 

The j)hugoid poles are now about half way between the origin and the closed 
loop short period poles (which at at s = -1.116 i- j 0.389 sec“^) as indi- 
cated on Figure VI-6. The damping factor of this location is 0.9 whereas 
the optimal selection provides only 0.7 damping factor. 

The thrust feedback dcj is determined as previously using the re- 
vised closed loop (on u/df) approximate dynamics matrix. Applying the 
same weight ratio for ^ of 0.0003 sec^/ft^ causes the h pole at the 
origin to move out to s = -0.179 sec~l without moving the oscillatory poles 
into short period territory, as indicated on Figure VI- 7. The resulting 
feedback gains are (in u,0, and h in consistent units of ft, sec, and rad): 

Ccj = (-.058, 3.895, .0173). 
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df 



FIGURE VI- 7 LOCATION OF CHARACTERISTIC POLES 
AFTER CLOSING FEEDBACK LOOP 
dcj 
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These results are tested, as before, in the five state model to observe 
the resulting true pole locations. The revised approximate solution is 
tested using OPTSYS to detemine if the mode interaction is decreased 
through judicious choice of pole placement by approximate modal analysis. 
The resulting true pole location is contrasted with the expected location 
on Figure VI-8, The apparent mode interaction is manifested in a 
decrease in damping of the oscillatory modes. However, no reversal of 
mode (pole) relations is indicated and the altitude (h) pole remains at 
the chosen position. The two modes remain separated, demonstrating the 
spectral separation assumption is not violated. The approximate solu- 
tion can provide a rough estimate of the closed loop response and a 
readily useable tool to determine feedback control gains without the aid 
of sophisticated computer programs; even when modes are separated by as 
little as a factor of two (as are the resulting closed loop poles in 
this case) . 

6,6 FURTHER APPROXIMATE ANALYSES 

The following discussion presents excursions into several violations 
of the separation assumption. The two mode (fast and slow) reduction is 
assumed to be valid while applying controls that are coupled to the fast 
mode onto the slow mode approximation. The first selection involves the 
elevator control. Second, the spoiler is applied to enhance the slow root 
configuration. In each case, the approximate solution is tested by apply^ 
ing its control laws to the five state model and evaluating the resulting 
closed loop characteristic poles. The analysis is carried out in the 
same manner as described in detail in the previous section; therefore, 
only the results are summarized in this discussion. 
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FIGURE VI-8 CLOSED LOOP POLE CONFIGURATION 
(REVISED APPROXItlATE SOLUTION) 
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The inclusion of the elevator as a phugoid damping mechanism is 
stimulated by its presence in the NASA work of Messrs, Franklin and 
Koening (Reference N-2), The elevator feedback included a pitch angle 
term in the stability augmentation used in that study. Although the 
application of elevator control to the approximate phugoid mode 
violates the concept of separation, since the elevator is included in 
the short period closed loop embodied within the approximate phugoid 
equations. 

The elevator is applied to the approximate phugoid mode to provide 
damping to the oscillatory poles. This is the role of the aft flap in 
the first control definition, leaving it free to be applied to the alti- 
tube transfer function. One attraction of this approach is the throttle 
(cj) control may not be required to arrest transient sink rates (assuming 
of course that pilot intervention would shortly include the rapid addition 
of throttle). The control definition includes both approaches, with and 
withoiit throttle (cj) included in the feedback control. 

The phugoid oscillations are damped utilizing pitch angle (<9) 
feedback to the elevator. The transfer function associated with this 
control is: 

^ (s) = .008 s (s + 0.199) (rad/deg), 

de s (s + 0.0919 s + 0.151) 

Applying root square locus techniques provides a root locus as sketched 
on Figure VI-9. Pole locations of s = -0.45 ± j 0.2 sec“l selected 
which produce the follovring feedback gains; 
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FIGURE VI- 9 ROOT SQUARE LOCUS 0(s) 

de 
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= 0,239 (deg sec/ft) 

C© = -95.526 (deg/rad). 

Dropping the Cu term moves the poles to s = -0.43 + j 0.35 sec"l changing 
the damping factor from 0.9 to 0.8 v?hile economizing on feedback terms. 

Two control approaches utilize the elevator damped phugoid to 

develop control of altitude; first with the aft flap (df) alone, and 

second with the thrust (dcj) included. The aft flap alone control is 

determined hy applying the control to the altitude state (e.g. h is)) . 

df 

The results of this approach are contrasted with the full state inter- 
action on Figure VI- 10. Application of both aft flap (df) and 
throttle (dcj) is taken in the same manner as done initially, the aft 
flap (df) is associated with speed (u) and the throttle (dcj) with the 
altitude state. These results are sho\<n in contrast with the full state 
poles on Figure VI- 11. The approximate solution poles are sho\<m as 
closed boxes and the resulting full state (true) poles are denoted by 
X*s on both figures. The results indicate that interaction is signifi- 
cant but not overwhelming. The approximation, though violated, provides 
an indication of the control response and in performing the root square 
locus a weighting ratio of desired state to control is obtained. The 
relation of these vreight ratios is an indication of the control influence 
that is being applied. 

Application of the spoiler as an altitude or sink rate control is 
investigated, utilizing the i ;ot result as a starting point. The spoiler 
control influence is investigated V7ith the non-spoiler dynamics equations 
Since these equations do not change significantly between the two equi- 
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libritim condition (see the equations in Section 5,3). The approximations 
of separation are further stretched since the spoiler is principally a 
w control as is evident from the full state equations. The spoiler is 
viewed as an adjunct to the previous loop closures of de, df , and dcj 
assisting the dcj in h control. The results are presented on Figure VI- 12, 
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FIGURE VI- 12 CLOSED LOOP POLES OF APPROXIMiVTE 
SOLUTION WITH SPOILER CONTROL 
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VII. FIVE STATE MODEL ANALYSIS 


7.1 AUTOMATED ANALYSIS USING THE FIVE STATE MODEL 


The control problem of determining state feedback control laws to 
arrest excessive sink rates due to engine failure is addressed using 
the full five state plant model and an automated control computing al- 
gorithm. The Stanford OPTSYS computer program computes the optimal 
regulator control solution subject to the system dynamic constraints 
by minimizing the weighted squares of selected states and controls. 

The program's capability is described in detail in Reference (B-1). 

OPTSYS computes the feedback control gains with their closed loop 
eigenvalues and eigenvectors. The closed loop dynamics matrix is 
provided as output also, making successive loop analysis possible with 
out the associated computational tasks. The feedback control includes 
all the states which can be a drawback if all the states aren't measur- 
ed. Therefore, OPTSYS control solutions usually require an estimator 
to determine the unmeasured states; v/hich the program can define (in 
steady state) given the appropriate input definition of measurements, 
measurement noise, and disturbance noise (see Reference B-1), Another 
approach (an engineering expediant as opposed to scholastic or scientific) 
to the problem of unmeasured states is to delete them. Perliaps the 
particular state feedback that is not measured is also not providing the 
major portion of the desired stabilization. If this is so, then that 
state's feedback, though required to provide an "optimal solution" is 
not necessary for acceptable control response. The scholastic approach 
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would be redefine the initial weights until the feedback in the un- 
desired channel was rendered negligible. However, this would be very 
time consuming since the current approach to minimize a state component 
in a particular eigenvalve mode is by trial and error. There is no 
readily available algorithm that can provide the desired eigenvalves with 
specified eigenvectors. 

7.2 THE ROLE OF EIGENVECTORS IN COST WEIGHT SELECTION 

Automated control solutions change the control problem from the 
computational task of choosing eigenvalves by determining the required 
coefficients in characteristic equations to selection (or guessing) the 
V7eights in the cost function that will provide the desired response. 

The presence of a ready computer and a control problem to be solved, 
often in a short time, provides a strong desire to obtai.n foolproof 
weights with negligible effort; v/hich on occasion may be a guess at some 
values. No foolproof weight selection exists to this author's knowledge. 
One rapid method to initially estimate the cost function x<reights is 
suggested by Bryson (Reference B-2). 

Bryson suggests that the designer identify the maximum expected 
value desired for each state of interest and each control. The tendency 
of the optimal solution is to spread the control so each part of the 
cost is treated equally in minimization. Therefore, weighting the squares 
of the states so they are normalized relative to their importance V7ill 
provide the desired result. The suggested weights are the inverse of the 
square of the expected value of each state. 
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Another method is demonstrated by the analysis of the previous 
section. That is first analyze by approximate methods us'; successive 
root locus to identify weight ratios then define consistent state and 
control weights for the automated solution which includes the ignored 
interactions. Both weight selection methods are applied to this pro- 
blem. 


Once an initial set of weights are determined the computer computes 
the control and closed loop response characteristics that may not be 
entirely satisfactory. If the solution is unsatisfactory the vzeights 
must be changed and a nev? solution computed. The question of which 
weights to choose and by how much must be determined by trial. In 
cases, such as the one considered here, where the number of states is 
large, the search for good weights can be a major task. The uncertainty 
of which weights to attack can be reduced by utilizing the eigenvectors-. 

The eigenvectors associated V7ith particular eigenvalues identify 
the relative’ amount of each state that is included in the mode of that 
eigenvalue. Any good text on state variable system theory shows the 
eigenvectors make up the modal matrix which, in turn is the transfor- 
mation of the system (in its state variable form) to modal form (Reference 
D-2, C-1): 

X = Qx* 

where: x* is the modal state vector. The dynamics matrix can be trans- 

formed into modal form (diagonal in the case of distinct eigenvalues) 
and solution determined by transition matrix of the form: 
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X (t) =^(t, to ) X,(o) . 


The modal form of the transition matrix solution is: 

(t) = x*(o) 

where 

= Q"' Q . 

Minor algebraic manipulation provides the state response in terms of . 
its modal components: 

x(t) = x(o) 

x(t) = qAq-1Q2^c(o) 
x(t) = qA.x*(o). 

The term _A-x*(o) is the modal response to the system. Q is the modal 
matrix made up of columns of the eigenvectors associated with the 
particular eigenvalues of the system. Each column provides the state 
component multiplier to the modal response which transforms the normali- 
zed characteristic response into state units. The eigenvectors provided 
by the OPTSYS program, therefore, identify the state make up of each 
mode (eigenvalue or characteristic root). 

Selecting a weight to change, once an unacceptable closed loop 
response is determined from an initial selection of weights, is aided by 
considering the resulting closed loop eigenvectors. Since the eigen- 
vectors show the state makeup of each mode, it is expected that the: 
largest coefficients represent the dominant states of that mode. 
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Increasing the state v?eight in the cost matrix associated with a dominant 
state in an undesirable mode should result in increasing the stability 
of that mode (moving its poles to the left in the s - plane). However, 
the desired result is not guaranteed. If several states are dominant, 
then it is possible that increasing a weight on one of those state re- 
sults, not in an increase in mode stability, but in a decrease of amount 
of tnat state presence in the mode! Such action indicates that the weights 
not only determine the characteristic location of closed loop poles, but 
also the relative state orientation of the modal responses. 

The present approach of increasing state weights is a trial and 
error scheme which may either increase the associated mode pole response 
stability or decrease the amount of state presence in the particular mode. 
If the dependance of eigenvector orientation upon the relative weight 
selection could be identified then a new design approach to control could 
be obtained. Often control problems are associated with specifying a 
fast response of a few states while allowing others to be much slov^er. 
Inevitably, it seems, the fast and slow states are combined in a common 
mode making the control specifications difficult if not impossible to 
meet. If feedback control could be defined so as to eliminate, for 
example, the fast state presence in a slow mode through a proper selection 
of cost vmights the control problem would, in cases such as is described, 
be simplified and possibly more efficient. 
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7,3 INITIAL GOST WEIGHT SELECTION 

The choice of cost function weights for application to the optimal 
control solution process is made in three different manners. The first 
is to apply directly the weight ratios from the initial approximate 
analysis (Section 6, 1-6, 2), The next choice involves reviewing the various 
approximate solution attempts and defining a set of average cost weights. 
The last method is the application of Bryson's rule (Section 7.2). 

The first example of cost function weights are chosen from the 
initial approximate analysis. The ratios of state cost weight to control 
cost weight arise from the root square locus gain required to provide the 
expected pole locations. In order to apply these ratios one set of costs 
are chosen arbitrarily since no unique solution exists to determine both 
state and control weights. The root locus analysis only provides one 
ratio pei: state explicitly analyzed, so there are twice as many weights 
to define. The control weights are chosen arbitrarily to be: 

B . = 25.0 
cj 

®de " ^df " (deg"^. 

The state weights are determined from the root square locus ratios as fol- 
lows : 

Aq - 3500. Aq = 3.5 (rad"2) 

®de 

Au = 14.25 Au = 0.01425 (sec^/ft^) 



Ah = 0.0003 Aj^ - 0.0075 . (ft-2). 
Bcj 
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The successive root square locus technique doesn't explicitly identify 
the weights of the states that are not chosen to be controlled explicit- 
ly; therefore the state weights of and A© are made zero. 

Several approximate solution analyses were made which provides an 
overall background of state to control weight ratios. These are reviewed 
with engineering intuition applied to obtain a set of state V7eights. No 
interaction between states is considered in the cost function so the 
chosen weights are diagonal terms in the cost function matrices A and B 
as listed below: 


Au = 0.001 

■ Bcj =0.1 

= 0.005 

Bde = 0.01 

Aq = 75.0 

Bdf =- 0.01 


A© = 500.0 
' Aj^ = 0.0005 

(Where the units of each term are the inverse of the state unit, to 
which they apply, squared.) 

Another set of cost function weights are determined by applying 
Bryson's rule. A reasonable choice of control and state excursions is 
estimated to define the independent state and control weights. For ex- 
ample, select a speed variation (Au) of 10 ft/sec as reasonable, then 
the weight associated to it is: 

A^i = JL == 1 = 0.01 (sec2/ft2) 

u^ 10^ 

The following table relates the reasonable variations and their result- 
ing cost weights: 
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STATE- 


WEIGHT 


WEIGHT 

CONTROL 

VARIATION 

VALUE 


A 

u 

u 

10 ft/see 

0.01 

sec2/ f t2 

A 

w 

w 

10 ft/sec 

0.01 

sec2/f t2 

A 

q 

q 

0, 1 rad/ sec 

100.0 

sec2/rad^ 

A© 

© 

0 . 1 rad 

100.0 

rad“^ 

\ 

h 

10 ft 

0.01 

ft-2 

B . 
cj 

dcj 

0.2 (approx 157o of 

25.0 



nominal thrust) 



^de 

de 

10 deg 

0.01 

deg"2 

^df 

df 

10 deg 

0,01 deg~2 


7.4 OPTBLA.L CONTROL SOLUTIONS 

The three sets of initial cost weights as they are defined are’ 
input with the five state inodel to OPTSYS to obtain the optimal control 
laws. The resulting feedback control laws are shown on Table VII- 1 
and the respective closed loop pole locations are presented on Figures 
VII-l,VII-2, and VII-3, The characteristic response as indicated for each 
weight set comparing the results to the approximate work. 

The first chosen weight set results in the pole configuration 
shown on Figure VII- 1. Also shown, are the expected pole locations de- 
rived from the approximate analysis (Section 6.3), Recalling the results 
of that solution, the feedback gains that were determined in that analysis 
had true pole locations totally different from their expected locations. 
Indeed they had no apparent relationship with the analysis beyond being 
stable (e.g. located in the left half s- plane) . Now, the weights from 
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TABLE VII- 1 


OPTIMAL CONTROL laws' 

i = Fji + Gu 
= Cjc 

In each case units are • consistently ft, sec, rad and deg. 
1) Weights directly from approximate analysis 



- 0.012 

0.0079 

- 0.0641 

1.292 

0.0068 

c = 

4.523 

- 3.244 

85.241 

- 581.32 

- 2.101 


3.262 

1.262 

- 45.091 

263.35 

1.385 

Weighti 

s determined by review of the 

approximate 

analyses 


- 0.058 

0.210 

8.479 

- 1.384 

0.063 

C = 

0.062 

- 0.178 

106.58 

- 205.56 

- 0.083 

i 

] 

0.125 

0.167 

- 56.25 

115.13 

0.057 

Weight 

s determined by Bryson's rule 




r 




— 


- 0.039 

0.024 

■ - 0.080 

4.128 

0.017 

C = 

1.073 

- 0.519 

109.51 

- 218.71 

- 0.483 


1.607 

0.272 

- 54.47 

109.77 

0.274 - 
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® EXPECTED POLE LOCATION 
BY APPROXn-IATE SOLUTION 


2.0 


X OPTIMAL SOLUTION 


X 


1.0 




• 3.0 




- 2.0 - 1.0 
X 


-h -1.0 


“h -2.0 


FIGURE VlI-1 CLOSED LOOP POLES DUE TO INITIAL 
APPROXIMATE SOLUTION COST WEIGHTS 


s -PLANE (sec- 1) 


i.o' 
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FIGURE VII- 2 CLOSED LOOP POLES DUE TO AVERAGE 

APPROXIMATE SOLUTIONS COST WEIGHTS 
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FIGURE VII- 3 CLOSED LOOP POLES DUE TO BRYSON'S 
RULE DERIVED COST WEIGHTS 
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that analysis provide reasonable pole constellations with the resemblance 
to the expected pole locations (Figure VII-1). The mode separation is 
about that which was expected but with changes in response frequency and 
damping. This demonstrates that the mode interaction that was ignored in 
the approximate analysis was the major factor in providing results totally 
unlike those that were expected. 

The next example utilizes the intuitively averaged cost weights 
obtained from the approximate analysis experience. The mtajor issue in 
the approximate work was mode interaction and the inability to avoid it. 
One might expect the resulting weights to be conservative, in that they 
ask for tight control. Indeed, such a result is apparent from the closed 
loop pole constellation on Figure VII-2.. The short period is driven far. 
out into the left half plane and the slow roots are pushed onto the real 
axis. One might conclude this is overpowering control in view of the 
mild movement of roots of the first case. Viewing the control laws, 
though, contradict such a conclusion (Table VII-1). The feedback matrices 
for the first and second cases have coefficients of the same order of 
magnitude. Indeed, most of the second case coefficients are. of lower 
magnitude than the first. The significant difference from this trend 
is in the thrust control channel, where the second case feedback gains 
are larger than the first case. The difference in the resulting pole 
constellation in this second case is due to a fortunate choice of weights 
causing tighter control without significantly larger feedback gains. 
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The application of Bryson's rule results in the closed loop control 
pictured on Figure VII-3. The pole configuration appears to lie somewhere 
between the high degree of control of the second weight set and the re- 
latively mild pole adjustments of the first weight wet. The advantage 
of this approach is it requires no analysis experience with the plant, 
one merely specifies the degree of each state that is allowed and also, 
each control, enters OPTSYS obtaining the appropriate control law. 
Occasionally though the results aren't as acceptable as are these, leav- 
ing the engineer with the issue of weight iterations described in Section 
7.2. 

Each of the optimal solutions have reasonable responses. The second 
weight set provides the most stable response which might indicate its 
control would be large. On the other hand, the first case with its mild 
solution v/ould be expected to require the least control application. 

The truth of this expectation is obvious from the control laws for each 
case (Table VII-1). 

The control laws show the three cases to be quite different than 
the relative pole configurations would lead one to expect. The relative 
gains in each channel are indicative of the control Output for a state 
disturbance input. The first case is most different from the last two, 
having the least amount of thrust feedback with much greater elevator 
and aft flap control. The last tv?o cases are much milder controls, In 
fact the last two cases are much the same with differences of an order 
of magnitude on, only, the speed (u) , altitude (h) and pitch rate (q) 
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channels. The largest difference between these last two examples is 
the feedback of pitch rate (q) to thrust. The second case has a great 
deal of this control, whereas the last one has practically nun“ (a 
difference of two orders of magnitude betv/een the two cases). The re- 
sult is surprising, as it is not apparent from the dynamical equations 
that the thrust would play a significant role in pitch rate control 
(i.e. pitch damping). . . - 

7.5 COMMENTS ON MODERN CONTROL SOLUTIONS 

Three different approaches provide three different control solutions. 
Each one optimal for the cost function applied. The application of modern 
control theory is demonstrated in these examples. Automated control 
analysis relieves the computational drudgery of hand solutions (which 
classical and approximate analyses tend to be) but does not relieve the 
control engineer. The focus of modern control theory is upon the definition 
or appropriate cost functions that provide controls that respond to speci- 
fications. The classical engineer must learn new terms and may approach • 
problems from a point of view that is initially divorced from the classical 
pole location root locus concepts. That is, until a solution is attempted 
without success; whereupon a correlation between cost function and the 
classical pole location must be developed to change the solution. Although 
this seems like a great deal of trouble, tt'jo advantages are derived from 
its use. First, complex systems are directly analyzed without eliminating 
significant interactions from the analysis. 
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These interactions may be beneficial for example: the pitch rate 

feedback to thrust of case two which provides the highly stable response 
without large control gains. Second, the control engineer can relate 
control to minimizing functions of state variables with the cost function 
architecture, and as such can cast the control problem in terms of system 
performance, a form which management understands. The approach also has 
one drawback, it defines control as full state feedback. The full state 
vector is not measurable in most real systems. However, modern control 
theory provides a solution, the state estimator. The engineer, not wish- 
ing to add complex compensation, such as an estimator, may investigate 
whether an acceptable control results from dropping the feedback terms 
of the unmeasured states in the control law. One may expect that this 
approach will work if the unmeasured states do not require dominant feed- 
back terms to provide the desired pole locations. This latter approach 
is investigated for tViis control study as it is attractive from an economy 
of parameters standpoint. 

7.6 ELIMIKATIHG feedback paths from OPTIMAL CONTROL 

The optimal solutions are evaluated with particular feedback gains 
eliminated to evaluate the need of an estimator. The approach is to 
first determine the optimal feedback gains v?ith the cost term on the un- 
measured state set to zero; then to drop the resulting feedback gains 
from the unmeasured states without changing the remaining gains and 
evaluating the resulting pole locations. The last two optimal solutions 
(case 2 and 3 of the previous section) are used. The states that are 
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readily measured are the angle of attack state (w) and the altitude state 
(h) . Altitude is commonly available as a measurement, on complex trans- 
port aircraft, however, it usually entails additional equipment beyond 
that needed for stability augmentation. The elimination of altitude 
feedback is attractive. The angle of attack, on the other hand, is not 
commonly provided as a measurement. Angle of attack metering systems 
are found on experimental and high performance aircraft with special 
needs. Inclusion of angle of attack (e.g. w) requires an estimator in 
the context of this study, so it is desirable to eliminate from the con-, 
trol law.- The angle of attack (w) feedback is addressed first. 

Before dropping the feedback gains; optimal control solutions are 
generated for the two cases of interest with the (w) weight set to zero. 
This is done to minimize the effect of making the w feedback gain zero. 
The characteristic response is then determined with the w feedback gains 
set to zero. 

The pole configurations, using the modified approximate weights, 
(case 2 of Section 7.3) show a dependence of the phugoid mode on angle 
of attack feedback (Figure VII-4) . The real phugoid poles (determined 
with the w weight included) become damped oscillatory poles with the v; 
cost reduced to zero. These poles move closer to the jw axis when 
the w feedback gains are set to zero. The short period and altitude 
poles are only slightly changed. The eigenvectors of these poles have 
a significant w contribution substantiating the hypothesis: 
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Phugoid Poles (A^ 0) 


Eigenvalue 


^1.645 


-0.350 


(sec"l) Eigenvector (units in ft, sec, and 

rad consistent with the 
states vector x) 

0.144 

-0.744 

-0.003 

- 0.002 

0.652 

-0.116 
0.276 
0.0001 
■ 0.0004 
-0.954 
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s- PLANE (sec-1) 
^ COST = 0 



--1.0 © A^^ = 0 AND 

w FEEDBACK = 


0 


-3.0 -2.0 



- 1.0 


1.0 



FIGURE VII-4 ELIMINATION OF x,; FEEDBACK 

FROM THE APPROXIMATE WEIGHTS 
SOLUTION 
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Phugoid Poles (A^ = 0) 

Eigenvalue (sec“l) Eigenvector Units in ft, sec, and rad 

consistent with the states 
vector jc 

0.603 ± j 0.517 0.088 ^ j 0.039 

-0.606 + j 0.503 
0.00007 t j 0.00005 
0.00002 t j 0.0001 
1.000 + j 0.000 

The cost function chosen for this case (case 2 of Section 7.3) causes 
the rqtation of the eigenvectors such that the traditionally separated 
mode states are mixed. The control law, in this case, suffers from the 
loss of angle of attack feedback causing a poorly damped oscillatory slow 
mode; that is not a traditional phugoid, but a heaving motion much like a 
ship at sea encounters. 

The response with the control laws determined with a cost function 
defined by Bryson's rule show a more traditional result. Removing the 
weight from the w state results in real short period mode poles and a slight 
increase in phugoid stability (Figure VII-5) . The subsequent setting of the 
w feedback gains to zero cause the short period poles to return to being 
oscillatory with less damping than they were initially. The phugoid poles 
move closer to the j(i) axis decreasing their* damping slightly. This choice 
of cost function weights maintains the standard aircraft modestate re- 
lations, and suffers very little in removing the angle of attack feed- 
back. The main result is a loss of short period damping from a factor 
of 0.94 to 0.7 damping factor and increase in frequency by a factor of 2. 
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FIGURE VII-5 ELIMINATION OF FEEDBACK CHANNELS 
FROM THE BRYSON'S RULE WEIGHTS 
SOLUTION 
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The phugoid damping is reduced from 0.67 t:o 0.55. 


This last case (case 3 Section 7,3) is investigated for sensitivity 
to eliminating the altitude feedback also. The altitude pole moves ex- 
pectedly to the origin and the phugoid moves toward the real axis increas- 
ing its damping factor to 0,80 (Figtire VII-5) . These results indicate 
that this control law can be implemented, without angle of attack feed- 
back, as a stability augmentation system (with speed control) which can 
be turned into an approach control sink rate alleviator by simply closing 
the altitude feedback loops when needed, 

7.7 SPOILERS INCLUDED IN CONTROL 

The spoiler control is a direct lift device which is apparent in 
its influence on the w state from the dynamical equations. The key to 
arresting sink rate growth is to be able to attain additional lift 
rapidly. The application of the only direct lift control may enhance 
the response to lift disturbances; such as the loss of thrust. In- 
clusions of this control is investigated to determine if it is effective 
in enhancing the closed loop response. The analysis is developed in the 
same manner as the three control investigation using the five state 
model and the OPTSYS optimal control computation. 

The dynamic model and initial weights are taken from the previous 
analysis. The spoilers are assumed to be deployed thirty degrees 
nominally during approach. The equilibrium flight conditions were com- 
puted for this case in Section 5.2. A comparison of the equations of 
motion about the equilibrium flight conditions with the spoilers de- 
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fleeted is not significantly different from their undeflected (nominal 
configuration) case (cf, dynamics matrices Section 5.)« 

The dynamics matrices of the nominal case are chosen to be used for 
the spoiler control analysis to determine if their effectiveness is 
Sufficient to warrant more careful investigation. Also, the initial 
cost function weights, as determined by Bryson's rule, for the nominal 
configuration are applied, as they appear to be a good compromise 
selection between plant response and control deflection. Additionally a 
small value control weight is defined for the spoiler to allow a full 
deflection of 30 deg (B^sp ~ .001 deg“2). 

The spoiler does not appear to be a significant control contri- 
bution according to these results. The pole locations, for each of the 
conditions chosen, are virtually the same as the three control case 
(Figure VII- 6) . The response is not changed and the control gaiiis are 
not significantly changed (Table VII-^ • The thrust feedback gains are 
reduced an average of 25 percent. The ether control gains are reduced 
somewhat also. The spoiler control inclusion takes some of the control 
burden without affecting the closed loop response. These results 
utilize the control in a linear fashion (the scope of this study). 
Perhaps a nonlinear control option is appropriate for spoiler 
control. 

7.8 NON-LINEAR SPOILER CONTROL 

The spoiler may be treated as a fast acting configuration device 
which may assist in alleviation of high sink rates. Plan a nominal 
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s-PLANE (sec-1) 




FIGURE VII-6 OPTIMAL CONTROL WITH SPOILERS 
INCLUDED USING BRYSON'S RULE 
WEIGHTS 
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TABLE VII- 2 


OPTIMAL CONTROL GAINS 
(Bryson's Rule Cost Weights) 
(Spoiler Control Included) 

(In consistent units of ft, sec, rad and deg) 

ALL STATES WEIGHTED 



-0.028 

0.021 

-0.041 

3.145 

0.013 

C = 

0.738 

-0.342 

108.27 

-190.07 

-0.386 


1.388 

0.335 

-54.91 

119.37 

0.308 


3.334 

-2.861 

13.83 

-433.88 

-1.792 

ANGLE OF ATTACK ( w ) 

SET TO ZERO 




-0.026 

0.018 

-0.102 

2.919 

. 0.013 

C = 

0.847 

-0.783 

98.07 

-209.08 

-0.394 


1.337 

0.375 

-53.31 

121.31 

0.298 


3.149 

-2.380 

21.96 

-404.03 

-1.789 
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approach with partial spoiler deployment which is immediately retracted 
upon the loss- of thrust. The equilibrium flight conditions foi* these 
two configurations indicates the strength of this concept (Figures V-1 
and V-2). The approach condition for a 6 deg flight path with 30 deg 
of spoiler deployment is approximately 83 kt at 50,400 lb of thrust. 
Estimating the sink rate at the equilibrium flight condition v?ith one 
engine failed results in a sink rate of 18.4 ft/sec with the spoiler 
retracted. The equilibrium sink rate at the nominal approach velocity 
is 19.7 ft/sec. The same comparison for 60 deg flap setting results 
in an estimated sink rate of 18.5 ft/sec (at about 90 kt) . The initial ' 
deployment of the spoiler is worth about 1 ft/sec sink rate. This system 
also requires an engine thrust monitor to signal the spoiler retraction 
and as such it v?ould not react to other disturbance inputs. The spoilers 
do not contribute a great deal to the direct control of sink rate. How- 
ever, there may be a synergistic interaction with the other control 
capabilities due to the increased speed of approach with the spoilers 
deflected, v?hich should be researched. 
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VIII. FOUR STATE MODEL ANALYSIS 


8.1 SINK RATE ALLEVIATION THROUGH INCREASED DAMPING 

The effectiveness on controlling sink rate vjithout altitude feed- 
back is investigated. The control concepts thus far have included altitude 
feedback providing integral control on sink rate and requiring altitude 
measurement or estimation. The elimination of altitude feedback is at- 
tractive if it can be accomplished without the need of angle of attack 
information since the resulting state feedback set as the standard 
stability augmentation measurement set of states. This approach was un- 
successfully researched by Messrs. Franklin and Koening in Reference N-2 
stimulating its inclusion in this work, to determine if optimal control 
theory can identify harmonious interactions within the plant. 

8.2 INITIAL COST LTSIGHT SELECTION 

The automated control computation (OPTSYS) is employed to the four 
state aircraft model to generate and analyze possible control laws. The 
four state model is generated by eliminating the altitude state from the 
nominal model. The initial weight selection is based upon Bryson's rule 
incorporating the state relationship to sink rate. The control set is 
the nominal three controls dcj, de, and df. 

The initial weight selection includes the objective of minimizing 
sink rate. The optimal solution determines control that minimizes the 


cost functions: 


Objective of the control is to minimize sink rate. The state weights 
are selected to reflect that objective. The sink rate perturbation in 
terras of the model states is given by; 

h = w + Uq a . 

The cost function weights the square: 

ii2 = w2 + 2 Uq w©+ Uq2© . 

Selecting weights in relation to allowable values of states according 
to this relationship may result in minimizing the sink -rate response to 
disturbances. The state cost matrix is of the form: 

Auu 000 
A = 0 Aw'? 0 Aw© 

0 0 Aqq 0 

0 A^w 0 A©0 

The weights Aw, Ax*;©, and A©© are chosen in relation to the sink 

rate equation, arbitrarily letting the state values be unity: 

Aww ■- 1 = 1. ft"2sec^ 

A©w = Aw© = 1 , = 0,0037 ft"l sec rad 

2Uo 

A0g>= 1 . = 0.0001 rad“2 

Recalling that the equilibrium velocity is approximately 135 ft/sec. The 
control weights are chosen to be the same as those of the five state case 
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3; the application of Bryson's rule case (Section 7,3): 


= 23.0 
^de = 

Bdf = 0.01 deg-2. 

8.3 OPTBylAI, CONTROL RESULTS 

The optimal control analysis is made using OPTSYS to obtain a well 
damped response and investigate the consequences of dropping the angle 
of attack feedback. The initial weights are applied to the optimal con- 
trol problem with the resulting response pictured on Figure VIII-1. The 
short period is well damped at a frequency twice the magnitude of the 
earlier results. The phugoid response is nearly undamped. The eigen- 
vectors associated with these poles identify the troublesome states: 


(Units Normalized 

Eigenvector in ft . Magnitude 
Eigenvalue (sec“l) Unnonrialized sec rad) '^'100. 


-4. 

88 i- 4.08 

0.071 


j 

0.0A6 

0.29 



1.000 

+ 

j 

0.000 

3.39 



-0.023 

± 

j 

0.035 

93.4 



-0.006 

dr 

j 

0.002 

2.9 

-0. 

164 + j 0.332 

1.000 


j 

0.000 

32.88 



0.007 

± 

j 

0.008 

0.36 



0.005 

"ir 

j 

0.002 

3.91 



0.002 

dr 

j 

0.014 

62.84 
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FIGURE VIII- I POLE LOCATIONS RESULTING 

FROM FOUR STATE OPTIMAL FEEDBACK 
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The eigenvector components are normalized by non-dimensionalizing the 
state components with Uq and "T then computing the component magnitude 
for unit length of the vector. Clearly the major components of the 
phugoid mode are u and ^ . The velocity state u does not appear in 
the h equation rendering it an independent phugoid state component of 
the cost function. Enlarging the cost weight on u by one order of mag- 
nitude induces more phugoid damping (Figure VIII-1). Previous experience 
with eliminating the angle of attack v/eight with the intention of eli- 
minating w feedback has not affected the phugoid poles (Section VII). 
However this case is different from those previous cases. This difference 
is apparent from the closed loop eigenvector state makeup. 

The enhanced phugoid mode control just determined has increased the 
interconnectivity between the modes. The eigenvectors are rotated in the 
states to increase the interaction of traditional slow and fast mode 
states. The eigenvectors are shown as before: 

Eigenvalue (sec"^) Eigenvector Units in ft, sec & rad 

consistent with the state 
vector [u, w, q, 

-4.86 + j 4.07 0.071 t j 0.046 

1.000 f j 0.000 
-0.023 -i j 0.035 
-0.006 t j 0.002 

-0.453 i j 0.265 1.000 ir j 0.000 

0.049 t j 0.023 
0.013 t j 0.005 
0.017 -Jr j 0.021 
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Compare the slow mode with the unnormalized eigenvectors of the previous 
case, the w state contribution is increased. One might conclude that 
the increased phugoid damping is obtained by increasing the state inter- 
action between the two modes. The next step in attempting to eliminate 
w feedback is to reduce the w weight to zero to reduce the feedback 
gains. 

The closed loop response of the system derived with the w weight 
set to zero is markedly different than the two previous cases. The 
short period roots move to the real axis and into the region of the all 
poles of the other solutions in this study (except of course the two 
initial cases of this section). The phugoid poles are significantly 
changed, too. Their frequency of response is increased with a slight 
loss of damping (0.68 damping factor). The optimal analysis the w state 
influence is dominant to the system response. 

The elimination of w feedback must have a significant impact upon 
the pole locations since the system is sensitive to w control. The last 
step of eliminating w cost weight was taken to reduce the w feedback 
sensitivity, which is apparent from the resulting control law: 



-0.003 

-0.000 

0.017 

-0.034 

c = 

-0.794 

-0.194 

52.41 

-68.86 


2.872 

0.206 

-39. 11 

86.40 

(In consistent units 

of ft. 

sec, rad 

and deg) 



The corresponding control lav? with the w cost weight included is; 


78 


-0.009 


0.075 


1.100 0.397 


C = 


-1.137 8.076 236.42 -33.61 

2.957 1.312 -56.31 94.39 


The w gains are reduced significantly, as are the dcj and de feedback 
gains in most channels. 

The closed loop response changes significantly with the elimina- 
tion of w feedback resulting in a well damped system and well mixed modes. 

The closed loop eigenvectors demonstrate that w is dominant in both 
modes : 

units in ft sec & rad 
consistent with the state 

Eigenvalue (sec“l) Eigenvector vector Qi,w,q,e]'^ 

-1.32 + j 1.29 0.030 ± j 0.006 

1.000 + j 0.000 
-0.007 T j 0.009 
-0.006 + j 0..001 

-0.541 t j 0.400 0.165 ± j -0.079 

1.000 t j 0.000 
0.000 + j 0.003 
-0.002 + j 0.003 


The change in pqle configuration from the previous case (with w feedback) 
is demonstrated on Figure VIII-2. The w seirsitivity is apparent, the 
phugoid poles move a large amount as the short period is destabilized. 
The resulting system appears to be well behaved with approximately 0.7 
damping factor on each mode. However as is apparent from the eigen- 
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FIGURE VIII- 2 ELIMINATION OF w FEEDBACK 

FROM OPTIMAL CONTROL SOLUTION 
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vectors state makeup, these are not the traditional aircraft modes. 

The vertical motion damping approach results in apparently well 
behaved systems. It is possible that these control laws may control 
the sink rate through damping. However the amount of damping required 
to eliminate major sink rate buildup is not yet determined. Testing 
these results by simulation is one way of evaluating their effectiveness 
with the additional feature of presenting control input histories for 
evaluation of their relative efficiencies. 
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IX. ENGINE MODEL 


The engine model is provided by input-output thrust responses to 
step throttle changes. The data are provided in the form of the analog 
output of the NASA simulation (Reference N-1) . A first order response 
model is defined for use with the control results in evaluating the 
thrust lag degradation of the system closed loop response. The raw 
form data are rescaled for presentation and future use. The thrust 
responses indicate that the engine is not accurately simulated as a 
first order system. Further work to incorporate these higher order 
responses is identified. 

The engine 'model source is taken from the NASA analogue simulation 
output for several throttle setp inputs at initial thrust levels of: 

507o, 70%, 80%, and ].00%, kn example of this output form is provided on 
Figure IX- 1, the time history responses overall of the data appear to be 
of high order (such as third order), but well damped with increased 
throttle command (idcj) . The engine failure problem is expected to 
require positive throttle commands, for the most part. The positive 
throttle direction is determined to be reasonably modelled by a first 
order resp»onse representing the thrust lag. The raw data form is used 
to determine an average first order model applicable to the thrust and 
throttle step ranges expected to alleviate sink rate growth due to loss 
of one engine (25 percent of the approach thrust) . The approach thrust 
setting used in the current NASA simulation is approximately 60 percent 


iNtz PAGE BLANK NOT fllMEm 


82 



FIGURE IX-1 TYPICAL THROTTLE RESPONSE DATA (RAW FORM) 


full thrust. Loss of one engine is a 15 percent total thrust level step 
decrease, which corresponds to a 5 percent thrust increase in the three 
remaining engines. 10 to 20 percent steps represent the range of ex- 
pected control commands. The following responses are read from the raw 
thrust response data: 


INITIAL THRUST 
(103 lb) 


FINAL THRUST 
(103 lb) 


TINE TO 95 % OF 
FINAL THRUST (sec) 


15.0 

17.5 

2.5 

15.0 

19.5 

2.5 

15.0 

21.0 

4,0 

15.0 

9.0 

3.5 

17.5 

19.5 

2.5 

17.5 

15.0 

0.5 

17.5 

22.0 

2.0 

17.5 

11.5 

1.5 

22.0 

17.5 . 

1.0 

11.0 

17.5 

2.0 

11.0 

15.0 

1.5 

11.0 

7.0 

1.5 

11.0 

6.0 

1.0 


The full step change ( 195 ^) is chosen to determine the time constant 
because the thrust histories due to a 10-20 percent positive step change 
appear to be nearly first order and the dominant control response is be- 
lieved to be characterized by an equivalent lag over these command steps. 
The above data are averaged considering three partions of the data. 
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The first since the estimate is taken using the 50 percent initial thrust 
data, nominal approach thrust level is near 50 percent (between 50 and 70 
in fact). The average of that response is determined; then the 707o ra.nge 
is included; and finally the full data set is averaged to determine how 
much difference exists (eg hov? special the subsets are). The time con- 
stants at Tg^y^ are determined for the first order system, 

I_(s) = a 

Tc s + a . 

The three average time constants (in sec) and a (in sec"l) are determined 
to be: 


'^low ■ 

tg^ — 1.5 

a Qs 2.1 

'^med • 

t95 = 2.3 

.a 1.3 

Tall • 

tg5 = 2.0 

a ^ 1.5 


The engine model is not a first order system as is apparent from a 
more close observation of the step histories. The initial approximation 
for positive thrust changes is approximately first order. However, the 
downward thrust steps are oscillatory indicating a second order component 
A closer look at the upward data indicate additional order characteris- 
tics. To facilitate more careful studies of the engine model, the raw 
thrust response data is replotted on an enlarged grid to scale up the 
response curves (Figures IX-2 to IX-5). 

The thrust history data is scaled up to make the unique response 
characteristics more apparent. An enlarged grid of the same structure 
of the raw data form is chosen. Data points are taken carefully from 
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easily identifiable points of the individual time histories and faired 
with the aid of a french curve. The resulting thrust time histories 
are presented on Figures IX- 2, IX- 3, XX-4 and IX- 5 for initial thrust 
levels of 50, 70, 80 and 100 percent* 

The engine thrust response is more complex than the initial first 
order approximation. The downward thrust steps are damped oscillatory 
responses. The upward steps demonstrate the influence of higher- order 
poles., also. One particular observation of the low thrust (50 percent) 
condition is an apparent time lag in responding to step inputs; which 
increases as the step coiamand is reduced (Figure IX- 2), The low thrust 
upward steps also demonstrate a slow first order pole, other than the 
initial time lag, influencing the last portion of the rise to final 
thrust level. These effects change to an apparently well damped 
oscillatory influence at the 70 percent and 80 percent thrust levels 
where upon an oscillatory overshoot appears to increase with increased step 
input (Figures IV-3 and IX- 4) . The nature of the downward steps vary with 

initial thrust level, too; demonstrating more damping at the 70 percent 
thrust level than at either end of the thrust range. The engine plant 
appears to be a high order system of at least third order with poles 
sensitive to the initial thrust level and the input step command. 

The influence of the engine model response variation with command 
inputs is an area of research not included in this work. The premise 
of designing the control without consideration of the unique engine 
thrust responses has been done with consideration of this model uncer- 



tainty. The dominant firist order response about the nominal approach 
conditions provides poles in the region of 2.1 to 1,5 on the s*-plane 
real axis. The control analysis is done with this in mind by attempt- 
ing to keep the closed loop response poles out of this region of in- 
fluence, and by weighting engine control heavily so the closed loop 
response is not heavily dependent upon engine response. The success 
of this approach is evaluated against the first order engine model by 
simulation analysis. 

The influence of additional slow engine response poles (not in- 
vestigated in this research) is an interesting topic for further studies. 
The effects of the upward thrust response poles is not expected to be 
dominant under nominal approach conditions or conditions requiring higher 
thrust levels. Although the 60 percent thrust data is not available, 
it is expected that the lagging response is a low thrust phenomena since 
the 70 percent and higher thrust conditions do not exhibit any lag. 
However, the amount of lagging influence must yet be determined at the 
nominal approach condition. This lag was not considered in establishing 
the first order estimate of the engine model. Recomputing the first 
order time constant more carefully using the re-plotted .thrust histories 
provides a value of 0.91 when averaging the rise times to 63 percent of 
the step command of the 50 and 70 percent initial thrust data (Figure 
IX-2 and IX-3) . The following values are determined to compute the 
average first order ' coefficient value: 
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INITIAL THRUST LEVEL 

50% 

70% 1 

STEP 
CHANGE 
OF TOTAL 
THRUST 
(percent) 

RISE 

TIME 

t63% 

(sec) 

FIRST 

ORDER 

COEFFICIENT 

a 

(sec”^) 

STEP 
CHANGE 
OF TOTAL 
THRUST 
(percent) 

RISE 

TIME 

t63% 

(sec) 

FIRST 

ORDER 

COEFFICIENT 

a 

(sec"l) 

+35 

0.6 

1.66 

+25 

0.6 

1.66 

+20 

1.7 

0.59 

+20 

0.8 

1.25 

-20 

1.5 

0.66 

+10 

1.0 

1.0 

-15 

1.3 

0.76 

-10 

1.75 

0.57 

-30 

1.75 

0.57 

-20 

1.5 

0.66 




-30 

1.7 

0.58 


The inclusion of the apparent lag slows the estiitiated engine response 
down into the airframe dynamic region on the s-plane. It is expected 
that the results would differ in view of the sensitivity of the air- 
frame poles. 

The control laws will generally command thrust decreases as well 
as increases in setting from a disturbance input. It is hoped (to be 
shown by simulation) the relaxation of control dependence upon thrust 
commands eliminates the overshoot of thrust command sufficiently to 
eliminate the oscillatory thrust response influence upon the system. 

It is possible that more efficient and fast responses could be obtained 
by including these effects into a control solution# Researching the 
control possibilities with a non-linear or more complex linear engine 
model including a detailed estimation of the engine's expected response 
is recommended for further development of a sink rate control system; 


92 




beginning with a third order identification study of the engine data 
and proceeding to control laws with an eight state dynamic model. 
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X. SIMULATION STUDY 


10.1 INTRODUCTION 

The premise of this control system analysis is tested by; simulation 
with a linear model. The assumption was made at the outset of this study 
to ignore the engine interaction with the aircraft plant in developing 
control laws (as explained in Section 1.1). The dominant engine res- 
ponse, under nominal conditions, lies near the plant closed loop poles; 
and, as indicated in Chapter, IX, may, in some cases, respond slightly 
slower than some of the fast short period closed loop responses. This 
assumption is tested by evaluating the performance of several of the control 
system candidates with the first order engine model developed in Section 
IX for this study. The simulation approach provides an additional result; 
in that the relative efficiency and actual sink rate performance is 
demonstrated so as to be compared with the s~plane analysis inferences. 

In this manner, the analyses performed to define the control laws are 
compared. 

10.2 LINEAR SIMULATION MODEL 

Linear model simulation is chosen to be adequate. The linear air- 
craft plant model is readily available as it is defined to perform the 
control law synthesis. The linear dynamic model is defined about the 
nominal approach path. Thex-efore, perturbations due to the response 
to an engine failure in the presence of a successful control law are ex- 
pected to be within the linear approximate constraints (approximately 
10 percent of the nominal conditions). Furthermore, the linear model is 



readily available; ince it was defined for the control synthesis analysis; 
and its transformation to a time si.mulation is straightforward. The 
linear approach does not require the introduction of additional sophis- 
tication integration schemes. 

The linear simulation model is produced by adjoining a first order 
engine model to the plant dynamics matrices and developing a state 
transition matrix solution for a fixted time step. The time step is 
chosen considering discrete information theory constraints upon sinusoidal 
time histories. The closed loop aircraft and engine dynamics have responses 
within approximately 5 rad/sec and for the most part within the engine 
response of 2 rad/sec (cf Chapters VI - IX). The selection of a time 
step of 0,1 sec results in 12 to 30 samples per cycle for frequencies 
of s to 2 rad/sec respectively. This is within the requirements of 
discrete systems theory (4-15 samples/cycle as discussed in reference P-1). 

The linear* transition matrices are defined from the state dynamics 
equations: 

X = Fx + Gj^u^ -f- G 2 U 2 

where: 

Uj = control input 

^2 ~ disturbance input (engine thrust) 

The engine failure dusturbance is modeled by superposition, applying a step 
input of negative thrust (cj) to simulate an engine failure. The plaxrt 
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nominal approach model is applied since the nominal path is to be main- 


tained and as such all deviations are perturbations. The transition 
matrix is defined by: 

2£n+l = :^2Sn + A ^In = A H2n 

where; 


§ 

T - Ft' 

Fi= f e dr 



T = time integration step (0,1 sec). 


The traiasition matrices are defined by series approximation (as describ- 
ed in Reference F-1): 


^ = I + FT + F^T^ + F^t3 + , , . , 

2! 3! 

p . = Fit + + f2t3 -f 1 ^i • 

L 2! 3! J 

An eleventh order series approximation is used. The simulation time 
histories are computed by a FORTRAU V computer program (a listing is 
included in Appendix D) , 


10,3 SELECTED CONTROL LAW CASES FOR SIMULATION 

The control laws that are selected for computer simulation include: 
the approximate analyses results; the full state optimal solutions; an 
optimal solution, with A^ == 0; and the four state control solutions 
(summarized on Table IX-1). The particular cases are numbered in associa 



TABLE X-1 


CONTROL LAW SUMMARY 

(Units are consistent ratios of the state vector x in ft, sec, and rad 
and the control vector ju in deg) 

Approximate Analyses (VI) 


(VI- 1) "Optimal" Case: 


C = 


0.175 0.000 0.000 13.762 0.013 

0.00.0 0.000 40.000 0.000 0.000 

3.186 0.000 0.000 105.260 0.000 


(VI-2) Revised Phugoid Damping Case: 


C = 


-0.058 0.000 0.000 3.895 0.0173 
0.000 0.000 40,000 0.000 0.000 
0.4946 0.000 0.000 119.900 0.000 


Optimal Solutions (5 State Model) (VII) 
(VII-1) Cost Weights Directly From Approximate Analysis: 


-0.012 0.0079 
4.523 -3.244 
3.262 1.262 


-0.0641 

85.241 

-45.091 


1.292 -0.0068 

-581.320 -2.100 

203.350 1.385 . 


(VII-2) Cost Weights Determined by Review of Approximate Analyses; 


C = 


-0.058 

0.062 

0.225 


0.210 8.479 -1.384 

-0.178 106.584 -205.563 


0.167 -50.252 115.130 

(VlI-3) Weights Determined by Bryson's Rule: 


0.063 

-0.083 

0.057 



-0.039 

0.024 

-0.080 

4.128 

0.017 

c = 

1.093 

-0.519 

109.510 

-218.710 

-0.483 


1.607 

0.272 

-54.470 

109.770 

0.274 

Bryson's 

Rule Weights With A^^^ = 0, 




-0.035 

0.022 

■ -0.142 

3.689 

0.016~ 

C = 

1. 212 

-0.979 

99.784 

-241.749 

-0.520 


1.501 

0.331 

-53.003 

114.765 

0.273 


Optimal Solutions (4 State Model) (VIII) 

(VIII-1) h Damping Cost Weights With = 0. 

’'-0.003 -0.000 0.017 -0.034’" 

-0.794 -0.194 52.410 -68.860 

2.872 0.206 -39.110 86.400 
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tion to their section reference on Table IX-1 and are referenced to in 
this section accordingly. The control laws use in the three standard 
controls dcj, de, and df. The spoiler control is not included; as 
the results of Section 7.7, indicate that the spoilers are not 
effective feedback controls in the nominal approach configuration. 

Each control law is discussed in the order of their development, begin- 
ning with the approximate analyses results, 

10,4 SIMULATION PRESENTATION F0RKA.T 

The significant simulation results are shown as time histories of 
the parameters of interest in a common format. This allows ready com- 
parison of one case to another. The parameters are provided in standard 
aircraft terminology as follows: 

AIRCRAFT CONTROL 

MOTION . INPUT 

sink rate - H + li (ft/sec) thrust coefficient - cj - T 

qS 

forward speed - Uq + u (ft/sec) elevator deflection - de^ H- de (deg) 

angle of attack - 5= w_ (deg) aft flap deflection - df (deg) 

Uo 

pitch, angle - O (deg) 

Two figures are provided for each case; the first includes the aircraft 
motion history; the second provides the control history. The scales are 
Identical where possible, to facilitate comparisons (scale changes are 
identified where they occur) . 
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The time history results are output every second providing a sum- 
mary of the motion. Fine oscillations may not be apparent in the re- 
sults but none are expected either. This data rate facilitates hand 
plottingjthe method employed in this work, 

10,5 APPROXIMATE ANALYSES CONTROL LAWS 

The approximate analyses defined two control laws of interest 
(Table IX-1), The initial analysis effort was an attempt to define an 
optimal control law using approximate dynamics equations. The result 
demonstrated that the special separation assumption had been violated,- 
However, the result, though not appropriately, predicted by the analysis, 
was stable with well damped slow mode poles. Since the validity was in 
doubt a second sub-optimal control law was defined by combining arbit- 
rary pole, placement with optimal (Root Square Locus) techniques that did 
not significantly violate the separation assumption. These two control 
laws are simulated, the results are shox-m on Figures X-1, X-2, X-3 and 
X-4. The two cases are remarkably similar. The first case (VI-1) shows 
the excellent damping that is identified on the s-plane pole plot (Figure 
X-1). The revised case (VI- 2) is slower than its predesseor (VI-1), 
but allows a slightly lower peak sink rate, 17.8 as compared with 18.2 
fps (Figure X-2). Their control r-esponses are similar also, where the 
slowness of the second case. (VI-2) is apparent from the thrust (cj) and 
aft flap histories (df) (Figures X-2 and X-4) , The control histories 
show that the sink rate is controlled mainly by the thrust input with 
an initial aft flap transient motion assisting the response. This re- 
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suit indicates the response in alleviating sink rate is most likely 
sensitive to engine response. The engine model, included in this simu- 
lation (Chapter IX), reaches 90 percent of a step command in 1.1 sec. 
Whereas the thrust command asks that the loss in thrust be recovered in 
about 4 - 7 sec, for these two cases (Figure X-2 and X-4) . The engine 
demands are not critical (which is consistent with the major premise 
of this Study) and may tolerate the band of engine response variance 
(shown in chapter IX). Also, the thrust history is, in the main, 
upward substantiating the engine model assumption of a first order 
model (Chapter IX). This last observation is true for all the cases 
presented in this section. 

The approximate optimal solution, although it violates one major 
assumption, provides a more rapid control. This case (VI-1) uses 
faster control imputs (than the revised case VI-2) to bring the sink 
rate -back to 90 percent of its initial condition in 6 sec; which is 
during a nominal altitude loss of 92 ft and incurring an additional 
loss of 22 ft. The sink rate increase is stopped in about 3-4 sec 
or in a vertical travel of about 50 ft. The optimal approach although 
not productive in identifying the response (cf Section 6.3) does maintain 
the nature of optimality by its slightly superior response without in- 
creased control inputs. These solutions, though approximate, have the 
advantage of specific feedback channels defined by the engineer result- 
ing in an economy of feedback paths which eliminate (disregard?) un- 
measured feedbacks that are not absolutely mandatory. There is no w 
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feedback, hov/ever h is required since it is controlled (Table X-1; 
cf Chapter VI ). 

10.6 OPTIMAL CONTROL USING THE FIVE STATE MODEL 

The control solutions of Section VII are evaluated with the first 
order engine model of Section IX. The first two cases are the results 
of applying the weights derived from the approximate analyses studies. 
The last four cases are derived from the application of Bryson's rule 
without benefit of analysis and as such typify design by modern optimal 
control theory (as discussed in Section VII). 

The first two cases, VII- 1 and VII- 2, are the result of applying, 
first, the particular cost weights that arise from defining the control 
case (VI- 1) and, second by defining average cost weights with the bene- 
fit of the intuitive understanding of all the approximate analyses. 

Case VII-1 is compared directly with its approximate equivalent (case 
VI-1). The state time histories are very different (Figure X-5 compared 
with Figure X-1). The optimal solution has taken the cost definition 
with the system state information together to produce a control which- 
epitomizes the apparent objective. The cost weights objective of the 
first analysis is to eliminate excessive thrust inputs, and the optimal 
solution does ! The thrust is not applied in any significant amount, 
using large amounts of elevator (de) and aft flap (df) (Figure X-6) . 

The result is large angles of attack (10 deg) and a likeiv^ise, large pitch 
orientation (5 deg nose upwards). These responses are beyond the 
limits of the simulation linear model as evidenced by the equilibrium 
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performance map of Section 5.2; the condition of 132 ft/sec (78 kt) and 
12 deg angle of attack has no equilibrium solution at the thrust level 
of one engine out (37800 lb). These results are not appropriate linear 
control solutions since their application violates the linear pertur- 
bation assumption of maintaining controls and states near their equi- 
librium values. However, the initial portion of the time histories 
remain valid, indicating a possible nonlinear approach to controlling 
these transients, since the sink rate builds to only 15.5 ft/ sec in 
2 sec With this saturating control. Immediate application of thrust 
may, in conjunction with this control law (acting in the presence of 
saturation), provide a nonlinear control that eliminates sink rate 
transients. Pursuance of this issue is beyond the scope of this 
study, instead the intuitively averaged cost function control law 
is offered . 

The time history response to an engine failure using the control , 
law of intuitive cost weights (Case VII-2) is sho\m to virtually elimi- 
nate sink rate transients on Figure X- 7 without large control inputs 
(Figure X-8). The time histories demonstrate a mild transient; a peak 
sink rate of 14,5 ft/sec, and recovery in about 3 sec (less than 50 ft 
of altitude is travelled during the transients). This control law res- 
ponds as the pilot should forcing the thrust back to its appropriate 
value in about 1.2 sec. These results are sensitive to engine response 
since the command asks essentially for full replacement of thrust as fast 
as the engine can provide it. . 


The results here in this case exemplify the premise that the engine 
response can be successfully ignored when it is slightly faster than the 
closed loop plant poles it affects. The pole plot (Figure VII-2) indicates 
poles on the real axis inside of 2 rad/sec from the origin, and short 
period poles at about 2,5 rad/sec. The engine response is at 2.1 rad/sec, 
providing the results on Figure X-7, indicating that the slow mode poles 
(although moved from their desctibed locations) do not interact detri- 
mentally to the system (i.e. the response does not become oscillatory or 
unacceptably slow). This case would be an interesting candidate for re- 
searching the effects of engine response and higher order modeling as 
described in Section IX. The next case considered, though, is the case 
of Bryson's rule optimal solution. 

The control laws developed from applications of Bryson's rule (Case 
VlI-3) are a compromise of sorts (the cru>: of optimal control is to com- 
promise in a least square sense) betv;een the two previous optimal cases 
(VII-1 and VII-2). The state histories sho\^7 a mild transient response 
to the engine failure (Figure X-9) like the approximate solutions case 
(VI-1 and Vl-2) applying thrust elevator and aft flap in quantities that 
appear to be an average of cases VII-1 and VII-2 control contributions 
(Figure X-IG). The significance of this result is that it is obtained 
directly by using the expected control response in the cost function 
weights ’(Section 7.3) , demonstrating that an efficient full state control 
can be defined from the designers performance requirements without explicit 
regard to the system characteristics. The optimal control solution utilizes 
the system charaGteristics in determining the specified solution. 
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The results of this case (VII-3) are better than the approximate 
case (VI-1) providing a peak sink rate of 16.5 ft/sec 2 sec after engine 
failure (Figure X-9). The recovery, however, is slower reaching 15 ft/ 
sec in 3 sec and recovering to 13 ft/sec in about 9 sec. This is ac- 
complished with 5 deg angle -of attack and a pitch up to less than 1 deg 
nose down. The control demands are significant using 6 deg of elevator 
and 4 deg of aft flap extension (Figure X-10), This control law is 
the basis for further investigations concerning eliminating the feed- 
back paths (discussed in Section 7.6), The cases leading to elimina- 
ting w feedback are discussed. 

The analysis of Section 7.6 proceeds through three steps of state 
feedback elimination to remove first w feedback and then h feedback. 

The objective is to investigate control laws that do not require additional 
states beyond those provided to the standard stability augmentation ‘ 
system (eg u, g, and 6 ). The basis for, these attempts is case VII-4 
(Table X-1) developed in Section 7.6 by setting cost weight to zero. 
The resulting control law is applied to obtain the time histories on 
Figure X-1 1 and X- 12. The results are similar to the earlier case but 
using slightly more control deflection (Figure X-12) to obtain slightly 
more angle of attack and pitch angle causing an improvement in sink 
rate response by recovering in 7 sec rather than 9 sec (Figure X-11); an 
improvement lost by removing the w feedback. The above control case is 
modified to eliminate w feedback, the results of which are shown on 
Figures X-l3 and X-14. These results are similar to the two previous 
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cases with responses much as predicted in Section 7.6. The sink rate 
response is slower reaching a peak of 18.5 £ps In 3 seconds and recover- 
ing to 90 percent in 8 seconds. These results are also much like the 
initial approximate cases in sink rate response but not in control 
response. 

The subop timal control case (with w feedback eliminated) provides a 
control that transitions to a new thrust trim condition whereas the 
approximate case maintains the initial condition. The attitude and 
control histories demonstrate this action (cf Figures X-1, X-2, X-13 and 
X-14). Case VII-4w (without w) trims up to 5 degrees angle of attack 
losing about 2 fps (Figure X-13). The initial approximate case (VI-1) 
drives the aircraft up to 3 degrees angle of attack in two seconds then 
returns to the initial attitude without changing speed (Figure X-1) . The 
controls are different, too. Case VII-4w uses the moving surface controls 
to obtain a trim state v/ith slov? engine commands (Figure X-14) . The 
approximate case on the other hand commands thrust relatively fast, 
supporting the transient alleviation with the movable surface (Figure 
X-2), The last case with its slow engine commands may be insensitive to 
slower engine response models, and may be the best solution to accommodate 
such cases. The loss of angle of attack feedback does not seriously 
degrade the control capability. Next, consider the loss of h feedback. 

The consideration of eliminating h feedback is simulated and the 
resulting state histories are shovm on Figures X-15 and X-16. The 
rationale for such a move is provided in Section 7.6. The objective to 
eliminate sink rate transieirts through the closed loop system damping 
provided by the control. It is apparent from the results that this 
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attempt at least is not successful. The equlibrium condition with one 
engine out at the nominal approach speed provide a sink rate of 18.8 
fps (Figure 5-3), The steady state result from this case is 19 fps 
sink rate (Figure X- 15) ! The possibility that the control can provide 
damping without measuring h is discussed in Section VIII, 

10,7 FOUR STATE VERTICAL DAMPING CONTROL 

The results of Section VIII are simulated. The time histories 
demonstrate that the control provided is insufficient to restrain the 
aircraft (Figures X-17 and X-19). The problem with each of these cases 
is the engine is not commanded to increase thrust (a consideration that 
was overlooked in design due to the attention to transient behavior con- 
cepts yet ignoring the regulator nature of OPTSYS solutions!) A response 
which is obvious from the control law (VIII-1 on Table X-1). An example 
of how the solution should look is provided. 


A four state damping control law was identified during the course 
of these simulation studies that successfully alleviated sink rate 
transients. The control law provides vertical damping with the inclusion 
of thrust response. The control law is similar to case VIII-1 (Table X-1) 
with larger feedback gains to thrust (cj) being the most significant 
difference, with corresponding decreases to the other channels: 


C = 


-0.0005 

-0.4420 

0.7303 


0.0537 

-0.4606 

0.0010 


0.0118 

46.4900 

-30.1850 


-0.0039 

-48.5600 

44.9900 


(where units are in ft, sec, rad, and deg). 
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The aircraft response to an engine failure using this control law shows 
that sufficient control of the four body motion states (without h feed- 
back) can arrest transient sink rates (Figure X-21). This task is ac- 
complished with minor increases in angle of attack and pitch angle main- 
taining a nearly constant speed. The sink rate changes from the nominal 
approach condition of 13,1 ft/sec to a steady condition of approximately 
15 ft/sec in 4 sec having passed through a peak of 16 ft/sec at 2 sec 
after engine failure (Figure X-21). This act is accomplished by command- 
ing the controls as shoxm on Figure X-22. The elevator and aft flap are 
commanded with the engine thrust to minimize the sink rate transient. 

10.8 SUM4ARY 

Several control laws are evaluated; their performance demonstrating 
the relationships of various states to the sink rate motion control pro- 
blem. The approximate solutions have a common characteristic -- no w 
feedback. Another case (VII-4) is also evaluated without w feedback. 
These cases although different in detail provide nearly the same sink 
rate response, one that is slower than all the other cases that arrested 
sink rate. It is obvious then that w feedback is required to eliminate 
sink rate transients, thus requiring an estimator control system which 
is somewhat more complex (and requires an additional design task for 
the estimator gains). This analysis indicates the sink rates can be 
held to 17 ft/sec peaks at 2 - 3 sec after engine failure and recover- 
ing in less 10 sec of flight (which corresponds to 22-42 and l4l ft of 
altitude respectively) without providing w feedback, but including in 
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every case altitude h state feedback. It is possible that these res- 
ponses can be made more rapid with further work of the nature discussed 
in Section 7.6. 

The optimal full state feedback cases provide the best responses 
demonstrating the capability to eliminate sink rate transients, but re- 
quiring state estimators or observers for implementation. The possibility 
exists that the excellent response ofcase VII-2 may be made poorer with 
the addition of an observer or estimator providing slower response. 

Such analysis is required to apply a full state control. The estimator 
work is left for future efforts, as the focus of this effort is upon 
control definition. 

The issue of eliminating the altitude feedback, the only externally 
referenced measure considered (eg fixed to another axes system than the 
aircraft) is demonstrated to be possible, but with w feedback. Attempts 
to remove the critical w feedback did not result in successful actions 
(i.e. they did not eliminate sink transients). The vertical damping 
control is a desirable control approach which requires an angle of 
attack estimate in order to be applied to the study aircraft. 



. XI. CONTROL LAW APPLICATION TO AIRCRAFT 
U,1 PERTURBATION CONTROL 

The control lawd developed in the course of this research are per- 
turbation controls requiring care for their implementation to the study 
aircraft. The differences are apparent from the follov/ing two control 
implementation schematics, the perturbation mode on Figure XI- 1, and 
the actual aircraft description on Figure XI-2, During the analytical 
analyses a perturbation model is employed which strips the nominal or 
steady state contributions from the parameters within the aircraft 
dynamics model.. The control laws apply to the error (perturbation) 
from the nominal or desired conditions. v-Then one seeks to use the 
resulting control laws, care must be taken to apply only the error portion ■ 
of each feedback channel to the control law gains. The simplest approach 
to this problem is to analyze dynamical equations which have only zero 
steady states, which is impossible with moving aircraft. Furthermore, 
this study maximizes the non-zero steady state terms by using a nominal . 
approach condition to define the dynamics equations providing Vq, Wq, 
and Hq states. Measuring devices provide the full values of the state 
parameters being measured. 

The recommended application of the control laws resulting from this 
research is diagrairaned on Figure XI-2. The diagram indicates the states 
that are provided by measurements already attributed to the study air- 
craft (see Chapterl) . Angle of attack is not provided. The "estimator" 
provides the function of defining the control variables for the control 
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law. The estimator block may contain algorithms ranging from: (a) sim- 

I 

ply subtracting the measured state data from a preset reference . to' pro- 
vide the perturbation state to the control law; to (b) full state esti-- 
raators with control law feedback and manual input control feed forward 
paths that may bypass the control law. 

The development of estimators for each candidate control law is 
not accomplished as part of this research. The complete task is left for 
future work (in addition to several more identified throughout the body 
of this report). However, one example is described to demonstrate the 
control application and bring out some of the issues involved. 

11.2 EXAMPLE ESTIMATOR DEVELOPMENT 

The example applies to the, four state damping case discussed in 
Section 10.7. The four state dynamics matrix is applied so the in- • 
puts are in perturbation form requiring that each channel be stripped 
of the reference values (e.g. subtracted) for input. This is not a 
significant drawback because, of the measured states u, q and 6 , only 
one has a significant steady component (e.g. u) . The q state is zero in 
steady state. The 0 state is small, -3.2 deg (0.056 rad) and may not 
severely hamper the system if made zero. Hov/ever, 0^ is described here 
to be subtracted from the measured state. The ramifications of leaving 
it (©<} ) zero are left for simulation studies not attempted here. 
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The estimator is a full steady state Kalman filter for which the 
OPTSYS computer program provides a stationary solution. The solution is 
developed and described in detail in Reference B-1. The form of the 
solution is as follows: 

X = Fx + u + G 2 w 

y = Hx + V 

X = X + K (;y - Hx^) 

u = G X 

where: 

X = the actual state vector- 
^ = estimated state vector 
X = an intermediate state vector 
- measurement vector 
v = measurement noise (gaussian) • 
w = disturbance vector with gaussian noise 
K = filter gain. 

The filter gain is described in terms of the error equation dynamics as 
a feedback term: 

X = (F - KH) X ; ^ 

The OPTSYS program solution utilizes the disturbance noise covar^ 
iance and the measurement noise covariance in the same manner as the state 
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cost weights (A) and control cost weights (B) to determine the feedback 
matrix K in the above error equation. The concern with estimation is 
that the noise covariances are define terms unlike the designers choice 
of A's and B's and as such are not freely varied. Also problems exist 
such as including insufficient measures or defining disturbances that 
do not excites the plants modes leading to no solutions. These issues 
are discussed at length in references B-1, B-2 and F-1. 

The disturbance covariance is treated in the same manner as the 
state cost matrix in Bryson's rule application to define an estimator 
(e.g. a diagonal matrix). Also, the measurement covariance matrix is 
described to include the measured variables u, q and & and values chosen 
that typify their possible measurement precision. The measurement pre- 
cision. The measurement covariance is chosen as the square of represen- 
tative one sigma variances in the state measure. The resulting choices 
are as follows: 


State 

.Disturbance 

Variance 

Measurement 

Measurement 

Variance 

u (ft/sec) 

100. 0 ft^/sec^ 

u (ft/sec) 

1.000 ft2/sec2 

w (ft/sec) 

100.0 ft2/sec2 

q (rad/sec) 

0.001 rad2/sec2 

q (rad/sec) 

0.01, rad^/sec^ 

(rad) 

0.001 rad2 

Q (rad) 

0.01 rad^ 




which provides the error equation response sho\'ai by the pole locations 
of the s -plane on Figure (XI- 3) , The associated eigenvectors are given 
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below; 


Units in ft, sec & rad 
consistent with the state 

Eigenvalue (sec-1) Eigenvector vector Qi, w, 

-3.129 dr j 0.50 0.017 t j 0.070 

1.000 f j 0.000 
. -0.016 t j 0.004 

0.001 t j 0.019 

-9.99 -0.997 

-0.080 

0.003 

- 0.001 ■ 


-0.741 -0.010 

- 1.000 

0.002 

-0.0004 

Most of the poles are beyond the response region of aircraft plant 
sufficiently far so as not to interfere, except the last one above. 

It is an angle of attack pole that lies in the slow mode region which 
will interact w^^ith the system. The effects of this interaction may be 
tolerable. Therefore, it is recommended that such a test be performed 
in the context of further estimator studies. The filter gain provided 
from this solution is: 





(Where the units are ft, sec and rad consistent with the associated 
state vector x. ) 

which may be implemented in a simulation. 




XIX. CONCLUSIONS 


Several control laws to alleviate the transient sink rates during 
an engine failure during approach are defined with varying degrees of 
control in accordance with varying degrees of complexity. Comparison 
of analysis by approximate mode decomposition and by the full state 
variable optimal control, provides understanding of the dynamics and 
an indication of the limitation of the approximate analysis. These 
state variable control analysis solutions result in successful alleviation 
of sink rate transients whereas’ the classical autopilot did not. The 
classical autopilot allows peak sink rates up to 22 ft/sec within 4 sec 
after an engine failure, reaching a steady state sink rate of 18 ft/sec 
in 8 sec (corresponding to approximately 130 ft of altitude travel). The 
state variable solutions demonstrate an average response of 17 ft/sec 
peak sink rate in 2-4 sec, recovering to the initial sink rate in less 
than 9 sec, with an attendant smaller altitude loss, on the order of 
120 ft. 

The more complex feedback laws perform better. The approximate 
solutions and suboptimal solutions (without angle of attack feedback) 
provide less sink rate alleviation. However, the full state feedback 
control laws are yet to be analyzed with the limitations of an estimator 
providing the state feedback, so their performance quotes are tentative. 

The full five state feedback cases using the average approximate 
solution derived weights and using Bryson's rule derived weights control 
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the sink rate transient to which a 15 ft/sec peak sink rate. The former 
solution provides the tightest control recovering in 4 sec after engine 
failure (about 80 ft of altitude) by depending heavily upon engine re- 
sponse. While the solution derived from application of Bryson’s rule 
used slower engine commands to recover in 9 sec (about 100 ft of altitude). 

The approximate solutions and the reduced feedback optimal solution 
(w feedback dropped) allow larger peak sink rates while recovering to the 
initial sink rate within 8 sec. A peak sink rate near 18 ft/sec occurs 
in each of these cases, which is less than that of the classical autopilot 
(22 ft/sec). These solutions are of particular interest because they use 
the same feedback channels as the classical autopilot with the one addition 
of altitude feedback. 

The possibility of providing sufficient control (within h feedback) 
in the aircraft to eliminate sink rate transients is demonstrated. How- 
ever, the successful solution does not develop from an. optimal solution, 
but from an ."intuitive gain changing happenstance result" solutipn! .The 
optimal analysis is not a cure-all. The successful four state case pro- 
vides good sink rate control Without altitude feedback, but unfortunately 
requiring angle of attack (eg w) information. This case builds up to 
16 ft/sec peak sink rate in 2 sec, recovering to a steady state 15 f t/sec 
sink rate in 4 sec (about 70 ft of altitude). 

Throughout this research topics beyond the scope of this work are 
described for future work that may improve the performance of an actual air 
craft control law. The results of this work demonstrate several successful 
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candidate control laws that may assist the pilot in arresting the engine 
failure sink rate transient. However, along' the way several extensions 
of this work are identified ranging from; the basic issue of the cost 
function weight selection relationship to the modal eigenvectors; to, 
the addition of more complex analyses to test the assumptions and models 
used in this work. The dominant suggestions that may provide better 
control are summarized here: 

o Eigenvector constrained optimal control solutions, 
o Nonlinear control studies including verification of linear 
controls. 

o Engine model extension and integrated linear control analysis, 
o Estimator development to provide unmeasured states and possible 
engine thrust identification. . 

o Integration of the longitudinal control with lateral control 
requirements. 

The optimal control approach to developing control laws for complex 
systems does not provide a direct association of the modal states to the 
co.st function. A first cut solution to a control problem may yield un- 
acceptable response results; where the desired eigenvalues are determined 
but where the desired fast and slow states are mixed in the modal eigen- 
vectors, Additionally, if one could specify the state orientation of the 
closed loop modal response, then partial state feedback solutions could 
be more readily defined, without the trial and error methods demonstrated 
in this Study. The need of an understanding of the relationship between 
the cost function, the closed loop poles, and their associated eigen- 


vectors is suggested as fruitful research to reduce the art content of 
solution iteration. 

The topic of nonlinear control schemes, although not included in 
this research, may provide a better solution. Two possible nonlinear 
control schemes for future investigation are identified in the course 
of this work; approach with spoilers; and saturating control. One possible 
control scheme is to fly the landing approach with the spoilers fully de- 
ployed, so they may be retracted upon failure. Linear control laws 
developed about nominally deployed spoilers did not shov7 significant 
differences from the nonspoiler cases. However, there may be a syner- 
gistic interaction with the other control capabilities due to increased 
nominal speed of approach with the spoilers deflected, which should be • 
researched. Control saturation is the deflection of the controls to their 
maximum extent, to arrest' the transients. In this case, deflecting the 
aft flap and elevators so as to immediately arrest the sink rate, while 
simultaneously commanding an increase in thrust from the remaining 
engines, to hold the desired steady sink rate. One such possibility is 
identified within the body of this report, and should be analyzed further 
to prove its worth. 

The engine; although represented in this work by its dominant 
positive thrust step response, as a first order system; appears to be 
more closely identified as a third order system with coefficients de- 
pendent upon the initial thrust, the size of thrust change command, and 
the direction of that command. The possibilities- of better sink rate 
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control utilizing a nonlinear or more complex linear engine model inte- 
grated with the aircraft motion states should be researched as further 
development of a sink rate control system. It is recommended that such 
work begin with a third order identification study of the engine data, 
and proceeding to control law development with an eight- state dynamic 
model. 

The optimal full state feedback cases provide the' best responses 
to eliminating sink rate transients, but requiring state estimators to 
provide the, unmeasured states for implementation. The possibility exists 
that their excellent responses may be constrained by slowly responding 
estimators. Additionally, estimators may provide an acceptable engine 
thrust failure identification, for nonlinear control application. 

. The scope of this study is limited to the longitudinal control 
problems associated with an engine failure. An engine failure requires 
asymetric laterial directional trim compensation that may detrimentally 
affect the sink rate problem. It is recognized that these aspects will 
be included in the future evaluation of an acceptable control solution. 
However, the possibilities of dominant interactions exhancing the control 
or estimator capabilities should also be researched. 



APPENDIX A 


DERIVATION OF THE PLANT EQUATION’S OF MOTION 

The following is a derivation of equations of motion used in the 
S.TiO.L. engine failure during approach study in state variable form: 

X = Fx + Gu 

The study is limited to longitudinal flight (e.g. Translations in x and 
z body axes and rotation about the y body axes,). 

From Newton's Law: 


F = ma = 


v?here 


V = 


U + u" 

= 


W + w 

II 

•N 

1 


differentiating 


V = 


u 

• 

w 

L 




u 4- qN 

w - qu 


where 


q = , pitch rate. 

The rotation is governed by 

M = lyy q* 

Expanding the forces and moments about a nominal approach condition in 
a Taylor's series provides the state variable form of the equation of 
motion: • . , 
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Reorganizing the last equations so the derivatives appear on the left 
of the equation, alone and removing the equilibrium conditions results 
on the following perturbation equations: 


r -1 
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dsp 


where the mass and moment of inertial terms are incorporated into the 
dimensional derivatives. 

The principal interest of, this- study is in the altitude loss due 
to major perturbations in engine thrust. The altitude relationship is 
defined in the study coordinates: 

h = W cos© + V sin© . 

Expanding about perturbations 

Hq + h = (Wq + w) cos ( 0Q -r Q ) + (Uq + u) sin ( 0^ + © ) 
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Hq + h = (Wo + w) (cos cos 0 - sin0Q sin0) 

+ (Uq + u) (sin cos 9 + cos sinS ) 

since ^ I 

Hq + h = (Wq + w) (cos - ^sin^p) + (Uq + u) (sin + 0cos 5^) 

The aircraft attitude rarely exceeds 10 degrees during approach so 
Qq ^ I is also true 

Ho + h = (Wq + w) (1 - + (Uq + u) (^ +0). 

The perturbation equation is: 

h = W + Uq ^ 

since, u * 9 q is second order small like 6 ‘ ©o* This equation is 
adjoined with the above equations of motion to form the equations of 
motion for this study. 
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APPENDIX B 


DIMENSIONAL AERODYNAHIC DERIVATIVES . ^ 

The dimensional derivatives aerodynamic contributions are 
developed. The aerodynamic coefficients include the effects of 
thrust on lift, drag and moments: 

cj s T 
qs 

qS 

qSc . 

oi, = tan'l j 

The forces are to be expressed as a Taylor series approximately , - 

(1st order) in the variables of motion about the trim condition: 

X = - D cosc< + L sin<^ 

Z = - L cos - D sin<=d 

M = My . 

First consider the prime aerodynamics (e.g. controls = 0): 

D = Cj5 (^, Cj , df) ^ f V2 S 

L = Cl (^, Cj, df) % pv2 S 

My = (©<:, Cj, df) % v2 Sc 
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where 

v2 = u2 + w2 

= (Uq + u)^ + w2 
= + 2UqU + + w2 

= Uq^ + 2UgU (to first order with Wq = 0) 

Since the thrust is included in aero coefficients j the appropriate 
variables of expansion are: 

u, w, Cj, q, and the controls (de, df, dsp) . 

Expanding the axial motion equation (the first equation); 

X = X^ + X^ u + X^lw + X 1 q 
o lo I o 

X = D cos 'x!.+Lsin<^ - 

The equilibrium conditions are: 

Xo = - cos -f sin 

Determine the first perturbation derivative: 

X- b cos + L sin ■) == -jp {cosc^'j- ^ ( cos 
^ u . <3 u ^ u 

+<?L (s incQ L £Csin«>() 

Bu du 


is typically small so 

differentiating the cosine: 



Therefore coso^ 
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sin/wj (i /w\= - w sinoi <^(u“^) 
[uj^inj ^ u 

= - w slnc^ (“l)u"2 

= wsin<>i=«>^sin<i^ 
u 

similarly the sine: 

^ sin<5C = cos^ w(-l)u"^ = - ^ cos << ~ ~ ^ cos »< 

3 u u 2 u 

The drag derivative is developed: 


_JL 

^ u 


(?)■ 


^ Cp % f v2 S 

du 3 u 

since 

v 2 = (u^ + V? 2 ) 

^ = Cjj^ % f v2 S + Cj3 ^ ^ S 2 U 
<J u 

The term is zero in subsonic flight, leaving: 


(^ = Cq e S U 

<? U 

The lift perivative is identically developed; resulting in: 


^ = Cl S U , 

Substituting the derivatives into provides: 
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For o(« 1. ^ ^ S U Cj 3 

The w axial perturbation is developed in the same manner: 

X = -^D cos e< - D cos ex? 

^ dvr d'w 

+3L s in oi + L s in i?< 
d\r ^ w 

The derivatives are determined: 


Id 

d w 


3D 
h W 


h COS /w y= - sin w 1. ~ " J: 

hot V u/ u u u uy 

= - JL sinc^ 

U 

sin = cos ^ =1 cos<^ 

<3 w 3" w U 

Substituting into Xw: 

X^ = - C e S W + % ^ v£ S Cj3 . cos - D sin 

U ^ ■ U 

+ Cl ^ S W + % S Cj^ sin + L _1 cos oi 

U ^ U 


d cos oi = 

) w 


= _A Cp % ^ V2 S = Cq % ^ v2 S + Cj) % ^ S 2 W 
(J w 


w 


Cp (0 S W + % e V2 S Cp 
U ^ 


= Cl ^ s w + Cr 1 % /o v*^ s 


V - ^ 

ix, L 
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«=ils: 




Including the definitions for D and L: 


= - (Cd 

^ s w + 

h 

u 


coso^ + 


S Cq sin«C 
TT 

+ (Cl 

^ s w + 

h 



sinoli' + 


S Cl cosfti 
iT 


The conditions of zero oi or imply = 

U 

yields: 

= - % f S Cd^ + % e v2 S ^ 
~ U 


^alsOj which 


S U Cd^ - Cl 


The lift equation is determined in the s^vf^e manner as above. 
The equations are: 

Z = - L cos o<. " D sin o< , 

differentiating: 


Zu = - ^ cos =<- - L cos c< - 5D sin c-t 
du du <5u 

2-w “ ■ iii GOSoC - L cos o<. - dP sintx, 
dw ' dw • dw 

Computing the derivates in the expansion 


= - Cl S p U cosoC - Cl q S c4 sin c<. 

U 

- Cp S ^ U sinfec + Cp q £ oc cos c< 

U 
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- sin c<. 
dw 

- P d sin ot . 
dw 





For approximation yields 

Z„ = - % p S U (Cl^ + Cj) 

The rotary force terms are negligible: 

X = 0 

q 

= 0. 

The moment equation, about the Y axis is: 

My — - Cj.| ^'Y.5 c j j h ^ ^ Sc. 

The expansion derivatives are: 

^^u = S 5 2 U = + ^ S c U 

M,, = + a Cm ^ p V2 s c + p S c ^2 

2) w w 

f S C + Cm % p s e 2 W 

^i<’U 

U 

Which for ^ % p U S c 
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The rotary damping derivative is significant: 


Mq = ^ 

>q 


= _l(%pV^ScCj^) 

dq 


= % p V2 S c 


> q 


The aerodynamic derivative Cm^ is defined as: 


Cmq “ • 

'2U/ 

= % P S c G„ c 


*q f ■ m, 


q 2U 


Allov/ing that the velocities are: 


■? 2 
V- 


Results in the rotary damping dimensional derivative 


M_ = % ^ U S c2 C, 




The control derivatives are straight fon'jard aerodynamic 
coefficient applications. The non-dimensional derivative coefficients 
are provided in the aerodynamics data and are dimens ionalized for the 
state equations as is sho\«?n for the thrust control ( cj): 


X, 


CJ 


^cj 


= % f S cos <?4 + sin«^) 


= % ^ S 1^- , cos oi. - ^Dqj 


in ad] 


Mcj = % f V2 S c C. 


m 


CJ 
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The derivatives of the elevator (de) aft flap (df) and spoiler (dsp) 

are determined by replacing the differential subscript (cj) with the 

appropriate variable. The aerodynamic coefficients are provided in 

APPENDIX C. Note that the coefficient C is taken to be zero, 

*“dsp 

since no data is provided for it in the aerodynamic data reference 



APPENDIX C 


The digitized aerodynamic coefficients are presented on the 
tables in this Appendix. The coefficients are referenced, uniformly, 
to the wing area (AREF = 1667 ft^) and the M.A.C. (CREF = 16.3 ft). 

The aerodynamic coefficients are presented for the range of angle of 
attack, alpha (ot) and thrust coefficient (cj) oyer which the original 
data is plotted. The angle of attack extends from -10 deg to +28 deg. 
The thrust coefficient range is from 0 to 3.85. Alpha is in degrees, 
all other coefficients are nondimens ional. The coefficients are 
arranged as follows; 

Table C-1 Aerodynamic Coefficients (Cl » ^D’ 

Table C-2 Elevator Control Derivatives (CL^eJ ^Dde» %de^ 

Table C-3 Aft Flap Control Derivatives (CL(if» Codf* %df) 


Table C-4 


Spoiler Control and Pitch Damping 
(^Ldspj ^Ddsp» %dspV^q^ 


Derivatives 



/ 


TABLE C-1 

AERODYNAMIC COEFFICIENTS 


CL 


CJT= 

0.0000 

. 0 0 

1 .££00 

£ . 9 0 0 1.1 

3 . £5 0 0 

- 1 0 .0 00 0 

-.10 0 0 

1 . £ 0 0 0 

3 . 0 C' 0 0 

3 .70'.itj 

4. £000 

-c.‘, OCh/v 

1.0 0 -0 .0 

. £ 0 0 0 

4.1z'00 

4 . 3 0 0 1.1 

-* • 1* U 0 *j 

•i.OOOu 

£■ . 1 0 Cl 0 

4 , 1 000 

5.3£i:i0 

Ti . £ 0 0 0 

5.3Ou'0 

. 00 Oil 

£ . 4 0 0 IJ 

4 . Z' 0 0 0 

5 .ac'Ou 

£• .30 00 

7 , C 0 0 0 

iC.OO'j'j 

£ . 5 Ci 0 0 

4 . zi 0 0 0 

6.1OU0 

7 . 3! 0 0 0 

S . 3000 

i-tiuOvO 

£ . '0 0 '5 0 

3 . 0 0 0 

. 5 OOO 

7 . 0 0 0 0 

5..' 000 

iS’ii . 00*0'.' 

1 . i 0 'J 0 

5 . £ 0 0 0 

4.50 0 0 

£• . 3 0 0 0 

>S . 1 0 0 0 


CB 





fiLPH.-! •:.JT= 

' j . '.1 0 0 0 

. £7 0 0 

1 .BiOO 

£.'.4 0uo 

3.3500 

-1 0 . 'OOij'O 

.£5 0 0 

- . 0 4 0 0 

-.50 0 0 

— .‘3000 

-1.1 00 0 

-i.O'jijO 

..oLi';>o 

.£500 

.1000 

- . 5 0 0 0 

- . 5000 

6 0 


. £ 0 0 0 

.£, coo 

. 5 0 0 0 

.•3000 

ic! . 0 0i."0 

. £• 0 0 0 

.3500 

1.1000 

I.I 0 0 0 

. 950 0 

1 =1 . 0 0 i.i 1.' 

. £5 ij i.i 

1 . £ 0 0 0 

1 .5 0 ij ij 

1 . 7 00o 

1.9500 


1 . 'I* 0 0 0 

1 . 3'jOO 

1 . 0 ‘J 0 

1 .950o 

£ . 4 0 0 Cl 

. 0 0 h U 

i . 1 0 0 0 

1 . 4 '.1 C. 0 

1 . ' 0 0 0 

£ . 05 0 0 

£ . 3 1 j 0 Ij 


fiLPHP ‘ ' J 1 = 

0 . 0 0 0 'j 

. 97 0 0 

1 .9;;uu 

£ . '3 0 0 IJ 

3 .35 0 0 

- 1 0 . ij 0 0 0 

.3.500 

.3300 

. .0900 

-.£3i.ii.i 

-.Zi5 00 

-£,0o0'j 

.1500 

-.04 0 0 

- . £4 00 

— .46 0 1.1 

- . 7 ,0 UO 

tO i. 0 0 Ij 0 

- ,,£500 

- . 4 0 0 0 

" • -» X '.1 

- , 7 0 IJ IJ 

- . 3 0 0 0 

1£ . OOOo - ■ 

- . 4 ij U 0 

-.5£0jj 

“ . r* r ' J 0 

-.341,11.' 

-1.1 000 

13.0000 

-.3o00 

-.5£00 

V*-r?0U 

- -.33 U 1.1 

— I .£000 

£4. OOOo 

- . -£ 5 0 0 

-.45C‘Ci 

" ■ t-E (j 

-.31UU 

-1 .£iJOO 

£3.‘jOO’.i 

-.•£4U0 

- . 3, £ 0 0 

“ . 5 0 G 0 

- . 7'JiJiJ 

“1,1 iJ 0 0 
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TABLE C-2 



ELEVATOR CONTROL DERIVATIVES 



CLD£ 


- :,JT= 

0 • 0 0 u 

. 37 0 0 

1 .35 0 0 

-lO.vO'.iij 

. 1 

. 0 Cl 1 7 

.0 U 1 7 

-'■r; i '.1 '.1 ij V 

.Ol-c.:! 

. 0 0 r:£* 

.0111 

6 - iJ U'."J 

: . 'M y3 

• Cl Ids 

.0l3d 

1 C . 'j t.l 

.1*055 

. 1 1 o 

.Old-3 

1 cl . I.l u V ' i 

.013? 

.0137 

.0C3t. 

. UUUU 

- . 0 ij 1 7 

- . 0 05 1 

. 0 0 1 7 

del . i,‘ 'J VU 

. Cf 0 1? 


-.■:n:i.34 




C£I£.E 



dLPHft- 


On -J 0 0 0 

.r'?00 

1 .SI'OV, 

-10. 

■V0»'*v‘ 

. 0 05 1 

-.0017 

- . 0 l?d'? 

~’iT- • 

6 ‘j C J 

.Oud'l 

. 0054 

. 0C'4 .S 

c* > 

Id. 

0 i’J ‘j 0 

.00 3 4 
^'00 id 

. 0041 
. 0 055 

. 0 C:5 0 

. Oil;- '=-1 

IS . 

UL*V.v 

• U V z- *4 

-O'lfSO 

.00'3£ 

• 

U U 1‘ 

.0054 

, 005! 

. 0)073 


V* 'J 

■ .U0d4 

.0044 

.0055 



C:iLi£ 



LPhTi— — O' JT-- 

0\ u ‘J 0 0 

- r 700 

1 - S 3 0 0 

-1 0. 0 0 0 0 

"*.050 0 

**-05£;D 

•-.05.10 

-£ . 0 lOOu 

7 0 

^nOSrO 

*" • 0 r*c* 0 

c-.. 0)0 Oil 
1 2 . 0 I.l i; •-> . 

- . 0 0 

- ^ C’)55iV 

"■-Or.’ 50 
”-0550 

".1050 
“ . 05 3 0 

1 £i . ) j 0 0 0 

- . 0 1 0 0 

”.0440 

--05c'0. 

d4 . lO'ji.iij 

* 'jij'50 

”•0.150 

- - 0 3 > 0 

dS . )j ij.O 0 

.lOOlO 

*".Oc-UO 

” . 0 H3 0: 


SO'.iij 

5) . 35 )■) 0 

oi..p5 : 

.1 )-) 1 d 0 

(.1 1 3 0 

. U l 37 

Olds, . 

.I.I1S7 

0133 

.1.1 103 

0 )3 34 

* 1.) 1 'd O 

OOC'S: 

. Oidi.‘i 

001 

.0 Ij 3 3 


3 00 I.l 

3 .. 35 0 0 

Ij0d4 

- . 0 014 

0051 

. 0 u 55 

0 oil 5 

. 0 055 


• V 0 r* 

OUSd 

• *jOS5 

0 0 D 

. 0 OS 5 

0055' 

. 0 075 


300 0 

, 3. 35 00 

0 zi 3 t.i 

— ■ U4- j ' (_i 


■ U r-t' 0 

lldO 

- . 1 05 0 

10'3u 

- . 1 05 0 



iJOy t» 

-.0730 

0531.1 

- . 07dO 


( 
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TABLE C-3 


AFT FLAP CONTROL DERVIATIVES 


,CL Dr- 



0 . 0 0 0 0 

. 3 ? 0 0 

1 .3.300 

£ . 9 0 0 0 

■ 3 . 35 0 0 

. -10. 0 0 0 0 

.01 0 c 

. 0 3 5 0 

.0300 

.0700 

. 08 '5 0 

-•2 . ij 0 0 0 

.0?U0. 

.0£50 

.0450 

. 05'5u 

. 03'5 0 

. 0 0 ij 0 

1 •- . >j 0 0 0 

.0100 
. 0 r 0 0 

.0£50 
. 0£4 0 

.0400 
. 0 35 0 

. 0500 
.0400 

.0-5550 
. 04 5 0 

IS.OuCnj 

£4.0.000 

. 0 £ 0 0 
:. 0 £oo 

,0£50 

.0300 

. 0 3 0 0 
. 0£ 0 0 

i 0 35 0 
.0£50 

. 0 4 0 0 
.0 300' 

c!o . U 0 *j y 

. 0« 3 0 0 

.0300 

.0300 

.0£5U 

. o£0 0 


CDD" 





RLPHB— r- O-JF-- 

0 , 0 0 0 1? 

.3700 

1 , 3 3 00 

£ ,‘;‘0OO 

■3 . '8-5 0 0 

— 1 Cl . 00'->0 
—£ . 0 0 Ov 
. 0 0 0 0 
ii'.oooy 

1 o . 0 0 0 0 

.01 00 
. 0 1 u 0 

X 0 1 0 0 

. 0 1 0 0 

. 0 1 0 0 

.0350 
.0350 
. 03 30 
. 0-£5>0 
. 0£50 

.0550 
. 0400 
.03-50 
. 04 0 0 
.0400 

.0550 
.o5oO 
. 05 00 • • 

. 0500 
. Or:- 0 0 

. 055 0* 
. Oo Cl 0 
. 085 0 
. 085 0 
. 0 ." '5 0 

•24.00U0 

. 0 1 0 0 

: .0350 

, 0 4 0 0 

• 055o 

. 0 7 0 0 

£3 . 0 0 00 

..010 0 

. 0£5 0 

.0350 • 

. 0550 

.0700 


OnDF 





— - i;..IT = 

0 . 0 0 0 0 

.'3700 

1 .'3 300 

£ . ’? 0 0 0 

3. 13500 

-10.0000 

.1000 

. 1 0 0 0 

. 1 0 0 0 

.1000 

‘ . 1 0 0 0 

-£.0000 
3. 0000 

.0400 

.0300 

.0-3 00 
. or 00 

.0300 
. 0-300 

■ .0300 
• 050o 

.O3o0 
. 09 0 0 

1 S'. 0 000 
13.0000 ' 

£ 4,0000 

.0500 
, 030 0 
. 0500 

. 0 3 0 i.'i 
.0400 
.04 00 

.0700 
.05 0 0 
.0400 

. 0:- 0 0 
. 04 0 O - 
.0£0u 

. 0':> 0 0 
,04'j*.t 
O.OOC'C' 

£•3.0000 

.0300 

.03- 0 0 

.0300 

.01 00 

.0100 


14 3 



TABLE C-4 


SPOILER CONTROL AND PITCH DAMPING DERIVATIVES 


vLDIr 


t J 1 = 

G . 0 ■’! il ■ j 

.. j' 7 0 0 

l.'?0-:0u 

3 . 3 0 G 0 

3. 3 50 

“ 1 0 . 0 0 0 

0 . 0 0 0 ''i 


- . 01 00 

- . 0 05 U 

O.fJtiii 

. 0 'j 0 u 

- . GO 4i'; 

- . G O'S^ 

— . 0 Go 7 

- . Og^o 

-.003 

i , tj 0 'j ij 

- . 0 Oo F 

- . 0075 

-.OOSi;- 

-.OUc'3 

— . 0-1 

! £ . Li.OOu 

- . 0 G ? 

- . 0 J Uc: 


- . OU5i 

^.!j05 

i . CiCtCiM 

-.00)7 

-.'■'075 

-.0133 

--.CllOr" 

— . U '.'s 

c “t- « U *J U *.* 

GOi? 

- . 0075 

-.0133 

-.G1m4 

-.015 

di'.OOOO 


-.ooi£; 

-.0 017 

- . Olb'o 

- . 034 


CI(D7P 


Hri 

U 

■ u C' 0 


• 3* *’ 0 0 

1 . 

95 00 

£ 

• 9 0 Li JJ 

3 

• 

3 0 L» 

ICl.OGC'O 

u 

IJ i’ *j *j 

0 

, 0 c c* c* 

0 . 

00 Liu 

0 

. 0 0 0 0 

0 


0 0 0 0* 

-3'.00oo 

0 

tj */ 0 0 


. n C 0 1; 

0 . 

0 C- 0 iJ 

0 

- 0 0 IJ L- 

0 


0 0 0 

z - . 0 0 0 

0 


- 

, u004 


0 009 

- 

. Ou U'-i 

- 


0 Ou-9 

13'. 0 000 . 

0 

C Ki *j C' 

- 

• 0 0> OS" 


001 ? 

- 

- OlU 5 

- 


uOl 7 

is.ooo'j 

0 

0 'j h 0 

- 

. ij'l'i'l 

- . 

0 0^ i! 

- 

. 001 V' 

- 


! j 0 1 C 

£‘4 . OiiOO 

0 

00 00 

-- 

. Gy.-; 1 


Ll »li4 5: 

- 

• 0 0^* r 

- 


uor ? 

H.3.00':0 

0 

0 000 

- 

. OOC 3 


OCfl? 

- 

,00t‘3 ■ 

- 


U L' 5- 9 


US'lI' i'T' 


Hf(- — C JT = 

0 . G O u »?• 

.9700 

1 a 9 5 GO 

iE;!.90u0 

2 a :• 5 V 0 

lO.OijL-0 

Ll a 1-1 i.t Li U 

G a 0 0 -0 0. 

0 a 0 0 0 U 

0 a 0 0 0 Ll 

C< a U 0 0 0 

-■3'. GO 00 

I'i » n rj r; f. 

G • 0 0 0 0 

0 . 0 f? 0 0 

O.Ouuu 

0 a Ll 0 'J 0 

€. . 0 0 0 0 

0.0000 

G a 0 0 IJ C» 

0 a Lit. 0 U 

0 a U 0 U 0 

0 a u 0 c 0 

li3.0GuG 

0 . 0 0 0 0 

0 a L' 0 0 0 

0 a 0 0 iJ 0 

0 a 0 U 0 0 

0 a 0 6 0 0 

13, 0 00 0 

. 0 . G 0 0 0 

0 a 0 G 1 j 0 

0 a 0 0 0 

0 a IJ U U iJ 

0 a U t j ij 0 

3 4 . 0 0 0 0 

0.0000 

0 a 0 0 0 0 

0 a 0 0 U 0 

0 a 0 IJ 0 l.l 

0 a U 0 0 L» 

£3 . 00 GO 

O.LiOO?j 

0 a 0 G C* 0 

0 a 0000 

0 a 0 Lf 'j Li 

0 a 0 ULrO 


cnQ 






3LPHri — 

— i3JT= 

0 a 

000 0: 


. 7 0 U 

1.330 0 

5: .’9U0.U 

3.3dOO 

-10. 

0 0 0 0 

-17, 

0 0 ‘G 0 

~ cLL.'* 

. 0 0 0 0 

- 5-1 . 0 0 0 0 

-56 a'i'Ouu 

-'57 . 50 0 0 

”“C a 

0 0 i.t u 

- 53 . 

0 0 Ij 0 

-37 

, 0 000 

- 5 T a 5 u G u 

-37-.50G0 

-37 .5o0 0 

6 . 

0 0 0 Gi 

:> *' a 

*] 0 0 'J 

-'43 

a 0 G ij 

-.54-'jOULi 

-5 5 .500'..- 

-5'5.0000 

13. 

0 0 u 0 

”CLl‘ a 

0 U L* 


a 0 L* 0 0 


-4 f .50 0 0 

— 1 7‘ . 5 0' 0 0 

13. 

0 1'.' 0 'll 

•*t 

0 0 0 0 

-13 

.0000 

“S‘5: a 5 l-0u 

-53. Go Go 

-41 .5000 

c4 . 

0000 

■ O a 

0*j00 

~tr:5 

a 0 U U 0 

-43.5000 

-■■1 lO , uHioo 

-59 a UUL'L- 

5! 8 a 

0 u 0 0 

-13, 

5 0 0 

~ 9 0 

. 000 0 

, -4 7.0 GO U 

-5o. gOoo 

-55. 500 0 


ORIGINAIi PAGE IS 
OF POOR QUALITY 


144 




APPENDIX D 


COMPUTER PROGRAM LISTINGS 


The computer programs developed for this work are included in 
this Appendix. They are as follows: 

Equil librium Performance Map ........ EQPERF 


Calculate Dimensional Derivatives ..... EQNMOT 
Linear Simulation Model SIM 


The subprograms required by each program are included followin 


each main program listed above 


♦ ♦ FOPHRF *‘t 


FILE I ^)4150^>i‘AiRPJ 


TI'-'E J. C 8 J 2 


C AERO 


SUBROUTlMfc, EQPERL - ’ — 

'blMfJMSlON C0lTmrf»'sr»Atr7) .CJTCb) tACU 

DliV.Ei'iSlOr'.. T(*>)»AL(1C1) 

) TABLE INPUT 

PRiNT/f • IMPIJT AERO BY ADO DAERC » ' ' 

rEAD(5»/) ((CC0EF(IVAR,1A,IC) »IApI , 7) f IC=1,5) f IVAR s 

READ(Nf/)fAT aA)»l Ai:t*7) . 

REAbC5f/)(C“jf(IC),IC“lf'5) 

PPlUT/f i nuiCK SEE AT SEA LEVEL DaTA A^'D 60 DEC FL, 
PRINT/ f* ILPUT^COOND RhO AREP WT’ 

REAOCSiT/) RHOtAREF*N 
PRINT /r ’ OF DSP* 

REAPCS f/) DF»D5P _ 

'lN=b. ■ • 

SHW= 1 , 

DNOs.66 

T 0 s’^l •tJiFf.Tvi 
TAi = 4 

PRINT /j.t I^PUT3 TMRUSTJ„ 

' R E A D ( S » /) ( T C f ) \ I =I . S ) ' 

PRINT /? * riPlIT 10 ALPHA* 

READ(H»/N (AL(I) ♦X=:i ♦ lu) ' L 

"DO 10 ITN)>Fs' 

PRiLT/f 

PRINT /♦ *_T’-ir!! 3T ='fTCn.1 „• 

print/ t* alpha, DE cl CD CN CAM* 

pRlf^T/»’ CJ V(KTS) V(FPS)* 

DO n I As 1 , i n 

AALsAUCIA) 

CALL AFR0(CnEFTAT,GJT,AALL5 

></ (ARFF^ACC 1 ,_) ) 

INsQ : 

TRsiNfl ■ 

CJ=T(IT)/COS*ARrFi 

call AEROCCOEFf AT,CJT»AALtCJtAC) 

Vs3QRTCABS(DS^2./RH0n 

CALL *ODOTcV>TMn ,nt’nT^DFQ« TQ,Tr ) 

'nc5 = bNDfYV!F;COS(GAN)/rQB 

DF=^ 1 . T ( AO ( 3 3 4- AC ( 9 ) ’F-DF t AC ( 1 2 ) *PSP ) / C AC C 6 3 ) 

C L = A C (1 ) 4‘ ( A C C a ) =1; S M W ) >1; [) E t A C ( 7 ) * D i- +7i C C 1 9 ) ^ D vS P 
CD = AC (2) + (ACfBmHW) YDE+ACC6)TDP +AC(in^DSPtDCD 
C MsAC C3 ) f ( AC ( 6 3 ) ^DE + AC ( 9 ) 4!Df' fAC (12) tQsP . 

CAHa i,.t . »CD/CL .. . . ___. 

feAT'LsATAF (C.AN) 

QBNsv>l^;CGSCOA?‘i)/{ARFrTCL) 

IFCABSCQBA-QBI.GT. ,01 .AND. IN. LT J, o 3 Ob = CBO 5 GOTO 

V=3GRT(ABS(08^2./RH03 3 

PRINT /jALCIA) ?0EtCL»CD»C‘l»GA:M>F'b7,295 

PRJNT/f * * tCDf V/1 .69^Vr!N 

PRINT/, (ACCI) .ACtI + 3) tACa + P3',AC(I+9) , 1 = 1,3) 
CONTINUE ■ ■ 

CQiMTINUE. . . . . . . _ 

RETURN 

END 


ltl 3 ) 


\ page 


'PS POOR QUAI-ii 


146 


♦♦ AEPO 


"" iaur 6^B05^AiRFJ ■■ TlMtrWVHr ' 



10 

11 

R 

13 

... 

lA 


IH 

10 

16 

11 

17 


16 


19 


?o 

R 

?1 

13 

?.?. 

C I A 





PG 


26 


J 


i pB 



20 


30 


S U 8 R q U T J, N F : A £ R (?{ C 0 E r > A T ♦ C J T » A A I y C J f A C ) * 

"D I MENS i ON' C 0E‘r ( 1 5 . 7 f 5 ) » A f( 7 ) • C J T ( b ) f AC ( 1 3 ) 

MC0EFsJ3 

NCs5 

NA=7 -■ * ^ 

AU=AAL 

C J I s CJ ^ • . 

“TrrA\l.UT.ATCn)ALl=ATa3JPRINT/f » ALP LOW TABUE’ 
IP(CJI.LT.CJT(U)CJI=Cjt(l)JPHirT/»» CJ LOh TABLE* 
IF(ALI.GT,AT(MA))ALI=AT(riA-n?PPihT/f* ALP_HIGH TABLE* _ 
IPUJI,GT;CJTCNC))ejI=CJT(NG-l) »PPlNt/» » CJ HIGH JABL^* 

DO 10 IA = J »^iA 

„ T FJ AL.I.EQ. .AT(!A) IGOTO U . 

TF(AU.LE,AT(TA)) IA::|A«UG0T0 U 

CONTINUE 

CONTINUE 

(50' R IC = j ,HC 

lEfCJI.EO.CJTacnGQTO K? • 

JI CCJ !• L E.CJTflCn I C^I c >.1 ;G 0 T O JL3_ __ 

CONTINUE 

CONTINUE 

IC Al. LOWER EDGE OF VARIABLE SPACE 
f<AtA=(ALI-AtCIAy)/'(AT(lA*-l)«ATtrA 
RATC = (CJI«CJTac))/CCJTCIC + n«CJTCIC)) 

DO ?.I) XCFsl.NCOeP 

" cl s'( C OE F (1 C F T I A +T'» I C C OFF ( fc F ? I A i fc ) ) * PA t A + C 0 E’f (1 C I A » I C ) 

C H = ( c; 0 E F (I c F ♦ I A 1 1 , 1 C + i) - C 0 E F ( I c F« 1 1 A ♦ X C P in =*" R A T A + C 0 E F a c r f I A , I C f 
ACaCF) = (CP-C.l)TRATCtCl _ _ 

CONTINUE ■ ....... ..... 

RETURN • . 


1 





pile:; 6W13054AXRPvI 03j2^ 


I 3MBf?.0UtlNL- T cViT » Ol'&Tf OHO»TO.f Tr-J ' 


2 

3 

H 

5 

6 

_ . . 7 _ 


Dr,M£i'M$IOi>| V/TCS) »AMT(5) »ATTC5) 

DATA (VT(n»I = l»S)/0.tl 1 U6Sf .167. ^!7 ♦223,29 » .'534.94 

^CAHTa).I = lTS)/1.0f1.0Ji 

&. i .02. 1 ,033f 1 .066 

S./( ATT ( X) »I = l .5) /l.Of 1.037, 1.C7.1. 115 1 1.227/ 

no 16 

fi 


IP(VTa).C7T.V)G0T0 1) 

9 

10 

CONIlNHiE 

)\i) 


I “G ^ 

u ■ 

n 

irU,UE.1)G0Tni2 

y? 


RAT = (V-VTa-l))/(VT(I)«VT(I-l)) 

13 


AWR=TAMTm-A’lT(T-»1 ) 1 *RAT + AWT(I-l) 



TGrc(ATT(X)«ATT(I- 1))TRAT+ATTCI-1 j 

IS 


nWnT=TMTnNO*AMRTT/(TO^TM*TC;R) 

16 


. HKTi'RM - , „ 

1 7 

12 

D^’DT=DMO^fT^I 

) B 


PETURM 

1 9 


P M n ■ 



EON MOT 


FU.Ei faal505 + AlRPj 


tlPKr oai 


__ subroutine: e qmmot . 

““lUHtNSION CdFF03,7»S)fAT(7)»CJT(b)tACTl3) 

COMMON / V A RMO T /F ( 1 0 1 10 ) » 0 ( 1 0 f J 0 ) t UO t v^O t THE TO f C J 0 1 ( 

I »DFOfDSPO»VO,GAMO 

AFrO tabu INPUT ' ' . 

PRlfsT /,f Ar .)0 AFRO TAfJLE FROM AI D tUF.MFNT AERO* 

RE A D iS » a (C 0 U LI ^ A Ri iL» l OUAglt? ) qC = l »5) »IVAR = 

RE%>CMn»V)7AT(fA) ♦IA=1 »7) 
rEadcst/) (CviTurn., ic=i 

WRm:(^j»/HATCIA) »IA=1»7) ^ __ 

’ WRITE(6»7)CCJT(IC),IC=J 
DO 12 IVARs:lfl3 

PR TLI J »* T.VAR =*tT.VA Rr * lA 12 3 4 5 6 7* 

DO J3 fU) 

E’RlNT/tlCt fCOFFaVAR»IA»IC) »IA=1 »7) 

13 CONThlUE 

la CONTIfiUF 
RAD-57.295 

OE=32_,_2 ^ 

PRINT /,* INPi.iT COND RHO AREf CREF W lYY* . 
READC5,/) RUO,AREF»CRF.F,N, AIV 

AMsW/GE . _ . . . , ; 

’ SHW-l , 

DM0=.a3 

’10 = 21,000. _ „ 

TN = A 

PRINT/, ’ DF DSP* 
read (5,/ 3 DF.DSP 

PRINT /, * INPUT FI.I0HT COND THFUST V GAM* 
RFAD(5t/> T»V,GAM 

G AM?G.AM^A D 

AAL^d , . 

Qd=. 5^RH0TV1:^2 
WS-K-+SIIvaVA’l3 _ _ 

WC=U^C0S(GAM3 

FIND EQUIt.IGRIUM CONDITXON 

IN S;0 

IJ = 0 

TGsTANCGAM) 

__PRXNT /,*^FiN[) EQUILIB C0ND3 t _ ' __ 

PRINT /t 

9 CJ=T/fQB*AREF) 

__P R J N T /, * VA R AA L - : 

PRINT /,AaI,CJ 

CALL AFROCr.n£F,AT,CJT,AAl.fCJ,AC) 

V = SQRTCGB*2/Rrlfj) „ 

7 ^’cALL MDGTI V',T,ONDT,DUO,TO, TNI • 

DCD=DMDT + VtCQS(GAM)/(QBtARE'F) 

__ _P^' 1 f - 1 • AXn_HAC (9) »Dr-*-AC ( 1 2 ) TLSF )_/ LAGC.P 
” CLl=ACO ) + (AC(0)^SHN)^DE+AC;(?>^DF4AC ( 1O3T05P 

CDl = AC C2) + (AC C5)4'3Hi-l3 Hcr)E + AC(8) =i^C‘F + AC (1 { j^DSPtDCD 
CALL AER0CC0EF,AT,CJTtAAk-2. »CJ, AC) 

V=SGRT (QB^2/RH0) ■ 

CALL 'MpOT(V,T,D''DTtPMO.T0,TN) 

_ ■PCD=DMDt»V*cn3 CGAM)./JLQBMR£Q^ ; 

DF 2 = * I ( AC c 5 ) lac C 9 ) ■’S^DF lAC ( 1 2 ) *C Sf ) / ( AC ( 6 ) TsHw ) 

CL2-AC C 1 ) + C AC C A) TSHU) TDF + AC C7 ) +AC (iO ) t-DSP 
CD2-AC C?3 4 C AC CS)TGtiV;) *pr. + AC (8)L^[ F + AC ( i U^D.OP + DCp 

149 



1 


) 


♦ ♦ EO(i!>‘OT 




rritj b/U3ur>*AiKrj 


pc.u^«cci- i-cuaj/ 2 . 

DCDA=(CDj-C');a)/ 2 . 

nAL= 0 , 

Tr(ABS(DCDA) .HT. .UOl) DAL = ( //S/ (nBrAREf) -CD U /DCOA 
irCAALfDA.L .LT. -lO, .OR. AAL + DAL .GT, 'd5,) 

S-DAL=DAL + .l 

Ii:CABS.tP.A L3-. • LT . .. . Q Oi. _ A D D . ... AiJ 3 £ PCI A GT^ ,.00,1 ) IM=lu 

CAUL AL;RO(COEr,AT,r.JT»AAU+DALf CJ»AC) 

VsSGRT(GEv<<^?VRf|f3) 

r.A.Ll. '-‘DOT.£V,TiD'’DT»n-MO,TO,TNJ 

DCO = OV-Ot'+'V*COSCGA'’)/CQB-+ARE:F) ' 

Of: --}.+• C AC ( 3 ) + AC ( 9 ) TDp + Ac. ( 1 1\) 4ir Sf )7 ( AC ( 6 ) 1=SHU ) 

... .CL. =AC C r ).+ CAC CJiO-?K3Hii.).=KDEJ-AC(.7LTLr tAG (i.Q.).LaSri.,_._ 

COsAC Cr2) + (5 ) *SHW) ♦DL + AC (8 .V^Df tAC C U ) *DSPfDCD 

WRITE ( 6 f/) » VAR AAl DAL DCLA DCDA Ii\ Tc CD/Cl CL CD» 
wRXTE( 6 ,/) AAL,D.ALfUCLA,DCDA,IN,Tn,CD/CL,CL.,CD 
WRITE:( 6 «/)< AC-3 ‘ 

DO n 1 = 3 » 3 

^R.I TE,C 6 r/):..AC.4l.)-» AC (I ♦ AC.ClJ-iiLi-AC7.Iih9> 

,U CQlMTI.^Ur 

. IFCA.aSCTG-CD/C.U ,GT. ,l)i)2 .AWD, IN .L^, 10) 

. I AAL = AAL.tP.AL*' IH = Inf l ? GdIO 9 

.TFCABSCDAL) ,LT, .003 .AND, IJ.LQ.O) 

&GANl=CD/CL‘DD3=n3?QGsl. = {GOTO R 

.. - IF £ I J .EG . UDODG-i aB-OBl )/ (C.P/.CLrX’AFa.) .Lqn.g(.TG^-GA*M ) » Dan 6 tUG ; 

R WF riE ( 6 » / ) ’ OBNEWs • ♦ OB 5 1 J = c ? I M = i- ? GCITOV 
IFCIJ.FG.F ,ADD. ABSCTL-Cn/CLI .Gj, ,CU^’) 

RGB-CQBi7GB)/2. JWRiTLCo,/) t OfhJLI- » »QS ? IJ»3 i 10 = 0 5 GO TO 9 
Tra-M.LT.jni -!RI TE C 6 » /I • SOLUTIf N CONVERGED » 

.&‘>L =AAL.+DAU 

- PRI7T /t : ^ 

C COMPUTE EOUILIB COMD IN BODY aXLS 
V=3GRT(QPJS./RHn) 

GA!!cATAN(CD/CL)- - 

PRIMT /t' EGUILIQ vBTAG AXFS VO GANG ALPD CJO TH DC 
PRINT / f V » G A * > ^ R A D I A .A L ♦ C J » T * U L 

.X--EI.ND. _3T.Aa. nPP TVG : 

PRINT /.» define F '-AT-HIX AMO G HaTRXX* 

C CALL. AERO ETC, . 

CALL, AFRD:(.CQEF,AT,CJT»AAU-5.»CJ»Aa . 

VsSQiR TTQB + P/RHO 3 

... . call! MbOTCV»Ta>MDT»DMO»TO,TN) 

DC 0 =D YD I tv *CJia.CllA-M.>/.C.,Qa7 AXEEI ^ .. . ' ' 

CLl=ACCi)h(:AC(4)*-3HN')'^DE + AC(7)YLf + .ACaO) YdsP 
CDl=AC (?,) f (AC(B) ^SHvO^OEfACCBj ^LFfAC C U)»D5Pi-DCD 
C;Mi=AC(3) + (ACC6)-^3h!;)YDEtACC9)^XFJACa2).^D3P . 

CALL AERD(C0Er»ATtCJT»AAL+5,fCJtAC) 

V=SQRT(GB;h?/RHa) ' 

.cALL„NO.G.Lcy.».i.,.x;,’aT.^nMai.xa....ti:i>._ 

DCD^7)NDT=i'V=^'COS(r7Ari)/(OB^AREri 
CLP = AC(3 3 + (AC (4)^-SH.0TDELAC(7) F + ACCli')Yr;SP 
. C.D 2 =ACC 2 ) + UCC.5T*SHN)TDC-fACC0)7LF-tAC(U)*D3PtDCD .. . ... 
CM 2 SAG r5) + (ACC A) YDE + AC(9)^r F + AC(l 2 3YD3P 

C DERIV PER RADIAN 

■ - CLA = RADY CCLS’^GL 3) / (10 . ) A.....-.- 

CDA-R^'iD-CCDR-COn/m, 

C A = R ‘j 0 ^ ( C ■ ' ? ’* C 1 ) / J. u , 

CALL AEROfCOEftAT ,CJT,AAL AT ) 


♦ ♦ EQNMoT itt KlUti bVl3u5l^AlRPJ ■ T r^t ; uB 

117 __V-SQR.T(RB^H/RH0) _ 

"il B ■ ■ C A U. '‘U)0 f ( V'Vf » f D-’O , TO r T )' 

119 DCD = Di*‘DTtV4-C05(GAM)/(QEJ^AlH>') 

\?i) CUsAC(n + f AC(4)4;5HU)^^DE + AC(:7)-*'r.FtACaO)+DS^ 

1?J CD,l=AC(2)-» f.Ar;(5)^SHVO tDE+AC(K)*DF + AC({n+BSP + C)C;D 

laa CMlsAC(3)i-(AC(6)?r3HW)=^<Dr.4ACC9)tC.F + ACa2)t05F^ 

123 C,AU AFROcCOFFt ATtCJT*AA LfCJf.5.AC) 

"'"l2/4"' Vs3URT(QG4-2/rH0) 

125 CALL f'DGT(VTrfDMDT?P'’.OfTO»TM) 

126 DCD=DMDT=4V*C0vBCGAM)/(0a-1^ARFF) ^ 

127 ' CU2=AC(i) + Uc(^)4\^HV04^DE+AC(7)4^rFtAC(.l63'*DSP 

12B CD2=AC(2) + (AC(5)4;.3Hi*!)4<DE, + AC(B)*(.F + AC(U)-4'D3P + DC:D 

t 29_ ^ CF*^2~AC ( 3 U ( a C ( 6 ) r S ri W 3 » D F, •*• A C ( 9 ) • * r F f A C ( 1 3 ) 4^ C » P 

" ""i30 ■ cuc:jscl?.^ci.i 

13.1 CDCJ-CDa-^-CDl 

132 CHCJ-CH2^CM1 . . L„ 

■ 1 733' ■“ CALL AFRO(COEF*AT,CJtf AAL 

13A V=3GRT(Qf5-t2/RH0) 

I 3 S C A t, L ^'ID 0 T ( V .„T jO MD Tj PMO ? T 0 , T N )__ ^ 

!36" DCO = LL’H)ft"v4:COS(GAH)/(GFr?ARiT'^ 

137 DE =^l .%CACC33+ACC^))*DF + AC(,l?.)4^l'SF)/(AC(654:sH!,-n 

136 • CL sACCj.) f (ACCin'4;sHM)*D[-:tAC(7)*C.F + ACC J 0)4'D3P 

139 ■ ^ ^ CD -3^C(H)4 CACCBl^SHWi^^DE + ACCSj^H-FtAC n U^DSP + DCD 

14(5 ■ CM =AC(3)4-(AC(6HsHL)1^DE + AC(9)4=DF + AC(t2)*OSP 

_ C OUTPUT AFRO IK_® fiO C0NJ:R0L„ PEER DEG 

14’2' "'print 7»‘'» CL CLA CLcJ CIDE CLDF CLDSP* 

143 PRINT/»CL«CLAtCLE:Jo (AC(4)4^aH/^) f AC(7) tAC CUO 

144 PRINT /f' CN CDA COCJ CDDE CDDF Cf:D3P« 

145 ' PRIf;r/» CD»Ci:}AfCDCJ.CACC5?'4SHA‘)4ACCA) ^ACCll) 

!4b PRINT /»’ C‘*;CHA C-'CJ CM[)E CMpp CFPDSP* 

147 __ ___ pRItsT /i CM ♦CMa «CMC J t ( AC (6) 4SHR3 * AC (9) ♦ AC (iPj 

I4>i PRINT /.’ C'-O 5^»»7 cM33 

149 c coHPurE stab dfpivs eoulib in body axes 

150 • CA=:C0S(AAL/'‘^AD) 

151. SAcESINCAAL/RAD) ■ 

152 U=V^CA 

W^V »SA: . 1__ N 

154 " ' " THET-GAU+RAQ-^AAL 

.155 ALsAALi/RAD 

1S6 C FORCE MATRIX 

157 fcT, tyj;CLTAREP45cRHO;ni4=SA«QD4:AU*CA/U)-CDtAREF4^(RHQ'}=U>CA+' 

158 aQ34^AL4=SA/'!J> 

159 FQ 5g)sCDtAREr4=rnB4;GA/U^Ril04?W*CA)>.CDA4^AREFYQB»-CA/ u 

l’61>’ ' ' R + CUAREFV(QBl^CA/U+PTI07hvi!aA)'K;i7A 

161 F (.1 t 33 r XVi 

162 F(2,: «>KCLTARFF-*(REFOtU4^CA + QB.?AL4^SA/U3«CDT.AREFY 

163 &(RH0’YUt3A“O8*AL4iCA/Ln 

164 f (2»2)j:CL-PAREP*(QB4SA/U-R|i04^U4-CA5'-CLA-<«AREFT-9i:j:^(;;3^/y 

__ 165 _ S>»CDYAREFH^ fDB4^CA/UfRHO»Vji:SA3"CQA 4^AhEF TQ84^S A/y ‘ 

_______ - _. ^ 

167 ! F(3i 1 )~CM%RHQ*ARFF4:CRLF4!U 

168 I _ E<3.a)-CMA*0RYAREF*CREF/U + CM4=Riie*AFiEF*CREFty 

169 I F(3»33=ACn33*QB4AREF4^CREF4-4^2/C2.>f;U) 

■ 170 C CONTROL matrix 

171 : Otj a^qb^arff _ „ 

f72 bBAC-QB-fAREP-CPEF 

173 G( I ? 1 ) -(CIC J'4^S5'i-CDCj4-'CA3^qBA 

174 . G(2H. ) = ('• } . tCLCJACA-CDCJ’NGAI^DBA . 

■ -. : 151 


EQNMQT 


FILt; 


/ 


TiiC* OH: 


\7^ 


G.ntJJrC.'-lCJ^QB.'VC : 

176 


G ( J, T ? 3 =: ( ( AC c -n t SH W 3 A- ( AC ( 5 ) *31- 

177 


G ( ? » ? 3 = ( • 1 . + ( AC C A ) *SHW,) '*f: A- ( AC ( 5 

17R 


Ga .23 = CAr.(6)*SHW)’*08AC 

17^ 


G ( 1 » 33 = ( AC ( 7 -1 tSA-AC < 8 ) tC A ) =*QB A 

180 


G ( ? » 3) C - J . 1: AC ( 7 ) >*0 A-AC (8) **SA 3 

181 


G C3 T 33 =AC, t9 3J1QB.4C. 

1 BP 


G (1 » A ) = ( AC ( 1 0 3 +3A»AC ( 1 1 3 *C A ) '+08 A 

1P3 


G C ? f A 3 = ( - 1 . t AC ( 1 0 ) *C A-AC ( 1 1 3 +SA ) 

ISO 


G(3»A)=ACCi2)+OF)AC 

18G 


DO U 1 I = ). t? 

1 8 ft 


DO 110 l *2 

187. 


F ( I » J 3 =ai ».J.) /A*l . 

188 

UO 

CONTINUE 

189 

Ut 

continue 

190 


DO J. .10 J = .t »3 . 

191 . 


FC5f J3=F(3»03/AIY 

19?. 


continue 

19.'^ 

. . 

, DO, 1 1 6-. -l-JJ L i.? 

19A 


no 1 17 J=i la 

1 9S 


G(Ii,J)=G(I»J3/A’' 

\9b 

1. .1 7 

CQMTINUE 

197 

1 16 

COhiTlNUE 

1 98 


no iia j=j»A 

1 99 


GCif J3 =G.(3.»J3 /Ai y .. . 

?0t) 

U3 

C0NTP8JE 

801 


1'.RITE(6»/) *. 3TATE3 ARE 1' U Q i 

?Q? 


^^RITE C6i/) L : F MATRIX’ 

80 3 


V.RITE(5»/3 (Fa,n*l2i,33 

?0A 


^••RlTE(6i/3 (FC?.,I) ,1 = 1, »3 3 

?05 


■.xRITECdi/) CF.(3tl3 »I=j 1,3), . . 

?06 


WRITE(6»/3 ’ 0 »'ATRIx DCJ DE DF 

?0 7 


NRITfcai/) CG(i >13*1-1 PP 

P()8 


WRITE(6»/3 (G.C2.I3 *Isl>a3 

?09 


:mRITEC 6>/3 (GC3>I) *I = lf/i3 

?i(; 

C NOW 

HAOE F G MATRICES FOR PASSING Dt 

?Jil 


..U.Q = U.,. .. : ... 

?1? 


V‘iO:2VV 

?1.3 


THFTOiTHET 

?1 A 


,CJO=CJ . , 

?1S 


DEO=OE 

?19 


DF0=6Q.+Dr 

817 


,.DSPD.;;a,aP - - - 

? 1 8 


VOrV 

?19 


GANOiGAM+RAD 

220 


RETUR9 . : . . 

?21 


EMD 


^:EXT PROG^^A* 



18 

QUM,^ 


152 


I 


* ♦ s I f' ♦ ♦ 

Hit 

1 

Sl'BPOUTltiL .5iy 

' ? 

DIi>?tN3lCN y (9) .M(0) turxi 

3 

nlf-*{.MSlCh P(6fo) .0(9.6) 

u 

/v;.RMnT/sr(io*ic 

s 

R.DFO.osPC.vn.nAfo 

o 

C If iTIALl^r P’PUTS 

7 

G0532«2 ... .... . - 

P 

RE:=;2090wOC)O. 

Q 

WA0=57.29S 

H* 

c STATrS n_W p_thft h_^ 

t 1 

C ZtKO ^•ATKICE5 

1? 

no 10 Irl,6 

1 3 _ 

on IS j.=i»6 

~ \a 

F (I.J)=0. 

IS 


i M 

G(I.J)=0'. 

1 7 

W(I ♦ J)=0. 

IP 

rci.j)so. 

19 

__ PCI . J) =0. - 

?0 

O(I»J)=0. 

P.\ 

V,Q(I . J)=0, 

?.? 

■ 15 CO.NTINnr. _ - 

P3 

X(I )=0. 

7’X 

U(I )=u. 

25 

, LD(n=e. 


10 cONTi'^'.it: 

?7 

C ALL INPUTS APf IN FELT DECS 

?P 

C INPUT 5 OF F G UO .10 THE TO C 

P9 

C ARh PY CO*'''‘'CG FRO.*’ FOUMPTOM 

30 

DO C 1=1*3 

3 ' 

DP.9 J=1»3 

3? 

r (I »U)sSMi.J) 

33 

G(I . J)=5Gf I , J) 

3^ 

9 CO.\TlNUF ..... 

3S 

8 CON'TIfiUf 

36' 

G ( 1 » a ) =SG ( 1 ♦ a ) 

37 

GC^iajssafataj- 

3B 

G ( 3 , a ) =v3G ( 3 ♦ fO 

39 

F ( t .4)=G0 

a 0 

Fca. 3)5-1 . -. 

<4 J 

F(5»2)=l . 


F(5»a)=l'0 

/i3 

C JEtiG IMF. .THRUST _n.I5TUrjEA;4CE -t\ 


C AMD 1 state FfOINF RESPONSE 

as 

DO n 1=1.6 

a6 

aCI.I) =0(1.1) 

a? 

F(I.6)=G(I.l) 

aP 

11 COiMTiMUr. 

a9_ „ 

DO 12 1=1 .6 _ 

SC- 

G(I ♦ 1 )sO. 

SI 

12 CONTINUE 

S? 

G(b. 1 )s2. 1 

S3 

F(c.6)=-2. .1 

sa 

DT = . 1 

. .. 55 

... c primtcut ".atricfs of the ff 

SB 

PRINT/.’ ' 

S7 

PRINT/, • r NATHIX. 

SP 

CALL Ppr(Fl 


t,ai3ob*Ai«pj 




r,rcs 


Dt-'O DFO D5PQ vO GA'^’O 


f’ODF.L 


153 


/ 


3 




rillt 6ai305><‘AlRPJ 


TU’E; u8t 


1 


suBmiTitiE-.: 

2 


dimension X 

3 


DIMENSION P 

4 


COMMON /VAR' 

S 


&,DFO*OSPOtV: 

6 

c 

IMITIAEIZE INPU 

7 


_......GQ.= 32.„.2 

B 


RE =2 09 04 000 

9 


RAD=S7.295 

1 G 

c 

STATES ,U. w 

1 J 

c 

ZERO matrices 

\? 


00 10 1=1,6 

) 3 


.DO. 

14 


F(It J?=0. 

IS 



16 


G(IfJ)=0. 

17 


W(I »J)=0. 

IB 


G(I,J)=0. 

A 9 .... 


C U.Ji.X=-a, — 

?P 


QCItJ.)=Q. 

21 


WQ(ItJ)=0, 

22 ■ 

15 

; continue 

25 


X(I)=0. 

?4 


UCl?=u. 


H 


25, 

H6 

?7 

25 

?<) 

30 

3.1 . 

32 

33 
3^ 

35 

36 
31. 
3B 
39 
a 0 

/>2 
, 4 3_ 
4 ^ 

45 

46 

47 
45 
.49. 

50 

51 


UDcn=a,- 

10 COiMTIMUE 

C AI..L I^:pUT3 ARF IN FEET DECS SECS 
C INPUTS OF F G UO. 4C^ THEID CJn DEO DF:,0 . DsSfTQ. V.D 
C ARE B.V CO‘-'MON FROM FQfiNGEOH 
DO 6 l:=l 3 
- DO .9 . J = l-i3- 


9 

8 


FCI»J3 =SFa »J3 
GCI ♦ J)=:3Ga , J) 

CONT INUE 

CONTIfinE 
GUr4)3SGU»4) 

^.,G.(.2.» .4.),= SJ1,C.2. ^ — 

G(3»4.)s$G(3?4y 
FU »4}:=G0 

F:(4,3}s-U.... 

F(5f2)F:U 

F(5»4)-U0 

C..j;NG.I,ME. tH.RliaL-ai.S,mEaA.Ma53iA.rRlX--. 

C AND I state FNGIME RESPONSE FODEL 
DO n 1 = 1*6 
,wcii nsGcr*n 
F(I»6)sG(I, n 
CONTINUE 

G(1»U=0. 

CONTINUE 


1 1 


J2 


S3 

F(6*6)=-2*.l 


S4 

DT=,t 


55 

C-, PRl-fiTQUT MA-TRICi 

r,3 DF-:TM-E. 

56 

PRiUT/f' ! 


S7 

print/,' F ' 

‘lATRiX • 

SB 

CALL PRMCF) 



...RROBLEM— . 


FILE; 


TI^Fr 03 J? 


S') 


PRIFT/» . 

fill 


PRIFT7,i 6 ‘IATRIX' 

61 


CALL PA9’(G) 

6? 


PRXFT /• 

hy> 

C F0F9 

’ PHI c;a?' lga''! matrices 

64 


CALL PHIGAMCE»G,W,Pf ntWQtDT) 

65 


PRTFT /. » PHI A'V'D GAM MATRICEc 

66 


CAlI.' PR!’GP) 

67 


PRliVT/» 

63 


CALL PRL(o) i 

6'J 


PRIFT /» 

70 


PRINT /»i engine DlSTRUBAliCE M> 

71 


DO 981=1.6 , , 

7?. 


‘WRITE! 6 tV) ‘Kl.n .WQd.O 

73 

93 . 

COWTiNUr 

74 

999 

COrJTlNUE 

75 


PRIM ?» 

76 


PRINT /t? XNPliT CONTRO.L CtATHSi 

77 

20 

PRINT/d j.JtCnEEi 

76 


RFADCSd) IM.CT 

79 


TF( J..F a. 0) GO TO 21 

S') 


CCI, J)=rT 

81 


goto 20 

6?. 

21, 

CONTI 'N.IE : 

B3 


PRINT/. L.C matrix' 

84 


CALL PRi'd) 

85 


PRINT/ j» initial CONDS TO Rl'N* 

86 


PRIAT/.' ALT MFOUT '’AX Tl'’E' 

87 


RE ARCS./) H-)»hEOUT.THAX 

86 


RO = t;, 

69 


ALT=Hn ... 

9 0 


TCUT=-1 ,7[>T 

9J 


DTOWTil , 

92 


T = 0, 

93 


ALP = WG/U04<R40 

q4 


H D 0 T = V* OF c 0 8 C T H E T Q / R A 0 ) + U 0^31 A ( 

95 


WRITE (6,/)* output' 

96 ' ' . 


>vRTTE('6t/) ' T RG 

97 


WRITE (6 d) ' TOTU ALP 

98 


WRITE (6./) ' U W R 

99 


WRITE (6i/)' DCJ PE DF 

1 0 0 

50 

CONTINUE 

FOl 


TE(T,LJ»TOUT1QOT0^.997 

lOH 


PRINT }f 

103 


WRITEfR.lOl) 'T.RO.ALT.HO 

104 

101 

format (SF 10,4) , . 

1 05 


WRlTE(6f 101 1 Un+X(l ) .ALP.THETQ 

1 06 


’wRITEC6t i'J.1 ) X 

107 


WRIT E f Of 101)1! 

106 


T0lJT = TOUTfDT01JT 

109 

997 

continue ^ 

1 1 0 


call VEQM(Uf X) _ . 

1 FI 


CALL VXMCUtC) 

1 1 2 


CALL r)VNCPfQ*WOf XflifUD) 

113 


.,T = rtDT : 

1 1 4 


H D Q T = F C 0 S C ' E T 0 / R A D ) + U 0 F S I M ( 

1 15 


A L P - C >< 0 F X ( 2 ) ) / ( U C. + X C U } F R A t ; 

116 


RO=ROf VrHcnscGANO/RAP) Fni 


Al T 
TI-FT 
TKF.T 
CSF^« 
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♦ » SIM piLtV ■ ■■■ ■ mtirGBS 

117 HQ=HQ«Vp+-aJ/iIfiAMny^^A.^ : -- 

~”flR AUT-HQ-X(S) 

UQ IF(HO.l.T.HKOUT ,AiJD. HD ( 1) .KQ, 0 , ) ID (V ) =- . 2b**C JO ; TOUT = T-DT " 

120 lP(T.GT,Tf AXIOnTO : : : 

121 IF(HQ,GT.503GOTO SO 

122 996 COi^iTI-gUF 

J .25 AHaTiiEg. CASF TYPE 1 !. ^ 

12A RFADCSi/) J-GAf'i 

i?3 upa)=o, 

126 .DO .996..I = i,».6._ ■ — - 

127 X(I)-0. 

128 U(I)=0. 

_1 29. 996 C Q N XX.MU.E: ^ — — 

no • IFaGAFuNE.OGOlO 999 

X'iS END 
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if:$ PHTGA^*' +♦ 


1 


5USR0!JTIME_^.PhIGAMT£! 

p 


MDNITnR 

> 

CP G 

PROM rx^GLi ^ AND Q A 

il 


Dir'lEiMSICN F(6,6)tG(6 

5 


D I { iF N 3 1 Gi\ D I ( 6 » 6 ) » 2 E 

6 


DO 10 1-1.6 

7 


DO 10 . J2l .6 

P, 


DI(I.T)=0, 

Q 


ZERa.J)=0. 

10 


D(ItJ)=0. 

1 1 


C (I»J)=0. . 


10 

COMTIN'JE 

3 J 


DO 11 I = ^ 

Ui 


D I C I . I) = 3 , 

IS 

11 

CQNT.tMUf: 

J6 


CALI. .EOrCPtDI) 

17 


0 13 0 (•' “ ,1 1 1 1 t w 1 

1 B 


IF{R.GT, 1 )GOTv) ).0 0 

1 1) 


CALL E:QKa£3lP.*P0-. 



CALL see (TEMP »T) 

?J 


call LQN(L»TE^iP) 

P? 

,1. 0 Q 

corti^iue , . : 

23 


TN = N 

?a 


SC-T/TN 

2S., 


CALL EQi'.(TE.MP,FO...... 

?h 


CAUL 3Ci'(TE’'P,3t;) 

27 


CALL XM (TEMP, PI 

2B 


CALL ADN(TEMP,Dn 

2« 


CALL EQ*'(P. TE'P j 

3(/ 

130 

CChTriTiE 

3 3 


CALL EQ!’(T£MP.»E:.3 

3? 


call XM(TF.MP.G) 

33 


CALL EQNCO.TEhPl 



CALI EOf(TE’'P.El . .• 



CAUL XMdUhP,.-') 

36 


CALL EQf'’(FQtTFMP) 

37 „ 


RETURN...... 

3B 


end. 

39 

C MATRIX OPERATIONS RRO'D 


UE| 64l305.'f:AJ,RpJ Tl'-’ 

P-? Q.» - 


•. Phi AMq t;Af, T 13 SA^-PUt DT 
63 fP’(6»P) *C. (hf6) »>’K6^6) fV-Q(6?6) 

C6?6) tt.(6»6) fD(6»6) »CC6>6J »TPMP(6»63 - 



ri 


♦ ♦ f)YM m.F:r*i><^l30S=^AlRPJ " ' “ TIt'E'roiUr/ 

j SUBROUTINE DYri(p»nfWQtXTU fUD) ' ^ : 

UCbifUDCb) 

’3 DIHENSIOM PC6fN),Q(6»6)»WQC6f6)tX(6)»TV(6)»TV2(6)tTENP(6»6) 

CAU- VEQf-’(TV,X) _ . 

5 ' CAUL VXr-'CTYfP) 

6 CALL EQf^CTEMPtQ) • 

__ 7 _ L A a^vLoti f ly 2 » U ) . 

B' ' CAUL VXf'’(fYP»fLf’p) ■ • 

.« CALL VADM(TV*TV2) 

JO call E0N(TE''P,WQ) ^ _ __ 

IJ call VEGM(fV2.UD) 

12 call VXP(TV?.»TEMp) 

1 3 CALL VA PM ( ly ^ ^ ■ 

Yu ^ CALL VF(3MCX»TV) 

15 return 

16 END ... 






adk ** 




el 


3 

U 

5 

h 

7 

B 


FlU:; t.^^l3(;5 + AlRPJ ’ Tlf^L; OBat^ 

_ SU8K0UT.1 UF . Ap' 1 ( A » f} )„ ^ ^ : 

D iND-iMS I Of"! A f b 1 6 ) » B ( f P ) . 

Msb 

. DO 10 lri ,M 
DO 10 J = 1 

30 A(lTJ)=Aa.J3+BCl»J) 

.... RETURM, ■ 

EMO 


♦ EQM 


t _ _ . SUBflQU.IIHE.EaL!.,CA»BJ 

2 DlMfNBir-M ACB. 6 ) »R( 6 * 6 ) 

3 ^’=b 

a DQ i 0 I - ! ‘-V • 

5 ' DO 1 0 Ji.t .M 

6 10 A(ifv.i 3 .=na»sn 

7 ; BEX,URiJ...._ 

B 

P c VECTOR - matrix ope rat I cm S' REO’O 


4=4= SCM 


t 

2. . b I M F N S i C A C 6 • b ) 

3 

U DO 10 T=1 

5* ‘ • ■ DO U) J™1 

6 10 A(I?J)aX*AC!,J) 

7 _ _ _ RETUR'XF 

' emd’ * ' ' 






?. A ( 6 » 6 ) » B C C3 1 6 ) » 0 ( 6 » 6 ) 

i = 6 

^4 .. DO i-0 J = 

5 DO iO I = 

6 DdfJ3=0, 

7 _„00_9' 0=1, M i 

B ‘ q 0(lf J5=r (i » Ji+Ad.'O^BCOfJ) 

q ,10 coMTiour 

10 DO U 

n DO 1.1 J=:} tM 

1?. 11 A(I»J)SD(IfJ) 

13 Return 

la" EMD 


fM 


_1 

?: 

3 

h 

s 

6 

_7 

B 

P 

10 
J 1 


SUB ROUT I [ v F TM(A) 

D I f»i£ N 3 1 C N .A ( 6 1 6 ) , D ( 6 f 6 j 

M = 6 

DO 10 X=1,M 
DO 10 J = 1 *M 

10 

_D0 = 

bb W JsltM 
l i ACI .J)=:Da ♦ J) 
rfturu 
END 


pri-t 


1_ 

S 

r. 

3 

a 

5 

6 
7 
B 


SUBROUT INF ?R:1(A1 

■" DINE NS ION Ad. 6 3 . 

DO 10 J=lr6 

WRITE (6. 1013 CACJ.I) t dlfb) 
101' ~ FOR 'TAT (6F .10,0 3 
10 CONTINUE • 

RFT U.R-M 

'end 




1 

?. 

7, 

a 

s 

b 

.7 


f b«i3uy+AiRPj 

SUB R.P I J 1 1 1 ‘% . . V A D ' 'j C A ♦ B.) 

AC6)»EUb) 

M = 6 

DO 10 J = t 
10 A(J)=A(J)+R(vn 
RCTUR^.' . 


"VF(3M' ♦!«< ■ ■■. 

j ■ SUBPOUTID^ VEOyCAjLB) 

^ - hem 5 I 0 N "K( 6 ) » B ( b ) 

t; f' = 6 

U DO 10 J-l»M 

5 10 A(J)=B(J) 

6 RETU9’\ 

7 „ _E!'iQ - 




3 

4 

5 

6 


S U 3 R 0 IJJ I NO V -5 CiJ..C k I .X.) 
Dli'’tN.5l0f'' A Co) 

M = 6 

DO 10 J = 1 
10 ACsn=X^A(J) 

RETURN! 

. EfvD 






VX"5 


1 

a 

3 

u 


7 

<i 

1:0 

11 

I? 


S U B K Q i J T I .V X MI A « B ) . . - — 

'D-TMEMSICN ACfe) »BC6 ?o) t0(6) 
M = 6 • 

DO 10 
D ( J ) - , 

DO D E = 1 

9 D CJ ) -Di iJ 0.+A.Cai *[iC J ifi) — 

10 cOMiriuc: 

no 11 jsif'i 

11 A(J)::DCJ). . . 

RETURM 

END 
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